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Abstract 

We consider an infinite chain of particles linearly coupled to their near- 
est neighbours and subject to an anharmonic local potential. The chain is 
assumed weakly inhomogeneous, i.e. coupling constants, particle masses and 
on-site potentials can have small variations along the chain. We look for small 
amplitude and time-periodic solutions, and in particular spatially localized 
ones (discrete breathers). The problem is reformulated as a nonautonomous 
recurrence in a space of time-periodic functions, where the dynamics is con- 
sidered along the discrete spatial coordinate. Generalizing to nonautonomous 
maps a centre manifold theorem previously obtained for infinite-dimensional 
autonomous maps [Jam03j , we show that small amplitude oscillations are de- 
termined by finite-dimensional nonautonomous mappings, whose dimension 
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depends on the solutions frequency. We consider the case of two-dimensional 
reduced mappings, which occurs for frequencies close to the edges of the 
phonon band (computed for the unperturbed homogeneous chain). For an 
homogeneous chain, the reduced map is autonomous and reversible, and bi- 
furcations of reversible homoclinics or heteroclinic solutions are found for ap- 
propriate parameter values. These orbits correspond respectively to discrete 
breathers for the infinite chain, or "dark" breathers superposed on a spatially 
extended standing wave. Breather existence is shown in some cases for any 
value of the coupling constant, which generalizes (for small amplitude solu- 
tions) an existence result obtained by MacKay and Aubry at small coupling 
[MA94| . For an inhomogeneous chain the study of the nonautonomous re- 
duced map is in general far more involved. Here this problem is considered 
when the chain presents a finite number of defects. For the principal part of 
the reduced recurrence, using the assumption of weak inhomogeneity, we show 
that homoclinics to exist when the image of the unstable manifold under a 
linear transformation (depending on the sequence) intersects the stable man- 
ifold. This provides a geometrical understanding of tangent bifurcations of 
discrete breathers commonly observed in classes of systems with impurities as 
defect strengths are varied. The case of a mass impurity is studied in detail, 
and our geometrical analysis is successfully compared with direct numerical 
simulations. 



1 Introduction 

It is now well established that many nonlinear networks of interacting particles sus- 
tain time-periodic and spatially localized oscillations commonly denoted as discrete 
breathers. In spatially periodic systems, breathers are also called intrinsically lo- 
calized modes [ST88] in distinction to Anderson modes of disordered linear systems 
[And58j . The properties of discrete breathers have been analyzed in an important 
number of numerical works (see the reviews [FWQSj IVMZ031 IDLHMS04j ) and their 
existence in periodic systems has been proved analytically in different contexts, see 
|MA94[ IAub98l [SM971 ES98l lAKKOll IPanOSl [Fk95l I.TamOB] and their references. In 



the context of numerical simulations or experiments discrete breathers often denote 
a larger class of spatially localized oscillations, such as metastable states, oscillations 
with a certain degree of periodicity, or even chaotic oscillations interacting with a 
noisy extended background |IKSF04l IGLCOSj . Nonlinear waves of this type are now 
actually detected in real materials [5504l ISwa99l ISES991 IEH021 IManOGj and also 
generated in artificial systems such as Josephson junction arrays, micromechanical 
cantilever arrays and coupled optical waveguides (see references in |CFK04] ). They 
are thought to play a role in various physical processes such as the formation of 
local fluctuational openings in the DNA molecule |PB89t Pey04| , which occurs in 
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particular during thermal denaturation experiments. 

Beyond spatially periodic systems, it is a fundamental and challenging problem 
to understand breather properties in nonlinear and inhomogeneous media, such as 
non-periodic or disordered crystals, amorphous sohds and biological macromolecules. 
For example the interplay between nonlinearity and disorder can provide an alter- 
native interpretation for slow relaxation processes in glasses |KA99-00] . In quasi- 
one-dimensional media, moving localized waves interacting with impurities |KBS94l 
EEEQSEPMM], extended defects |TP96j or local bends of the lattice (see |(::Kn4j 
and its references) can remain trapped and release vibrational energy at specific 
sites. 

The modelling of thermal denaturation of DNA and the analysis of its local fluc- 
tuational openings, also known as denaturation bubbles, represents another problem 
where heterogeneity is important. In order to describe these phenomena, a nonlinear 
model at the scale of the DNA base pair has been introduced by Peyrard and Bishop 
|PB89j and further improved by Dauxois et al |DPB93] . The model describes the 
stretching of the H-bonds between two bases, in the nth base pair along a 

DNA molecule (a large value of x„ corresponding to a local opening). Each bond 
fluctuates in an effective anharmonic potential V and interacts with its nearest- 
neighbours. The model is described by an Hamiltonian system, and can be coupled 
with a thermostat to study the effect of thermal noise in denaturation experiments. 
This model accurately describes the thermal denaturation of real DNA segments 
provided their heterogeneity is taken into account |CG98j (for example, the dissoci- 
ation energy of AT and GC base pairs are different). The Hamiltonian of the system 
reads 

^ = E f + ^-(^") + 2 + pe-^(--+-)) {x^^, - x^r, (1) 

n=— oo 

where Vs„{x) = -Ds„(l — e~°"'"^)'^ is a Morse potential depending on the base pairs 
sequence s„ G {AT,GC}. The case p = yields a particular case of Klein- Gordon 
lattice, i.e. the model consists in a chain of anharmonic oscillators with harmonic 
nearest-neighbours coupling. For parameters corresponding to real DNA sequences, 
Langevin molecular dynamics of ([T]) have shown that some locations of discrete 
breathers heavily depend on the sequence and seem to coincide with functional sites 
in DNA |KRBCU04] . but at the present time this conclusion remains controversial 

|vE(:;LPn6j . 

From a mathematical point of view, Albanese and Frohlich have proved the exis- 
tence of breathers for a class of random Hamiltonian systems describing an infinite 
array of coupled anharmonic oscillators [AF91j (see also the earlier work [FSW86j of 
Frohlich et al concerning quasiperiodic localized oscillations). These breather fam- 
ilies can be parametrized by the solutions frequencies, which belong to fat Cantor 
sets (i.e. with nonzero Lebesgue measure) of asymptotically full relative measure 
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in the limit of zero amplitude. These solutions are nonlinear "continuations" of a 
given Anderson mode from the limit of zero amplitude, and the gaps in their fre- 
quency values originate from a dense set of resonances present in the system. For 
disordered Klein-Gordon lattices, complementary numerical results on the continu- 
ation of breathers with respect to frequency or the transition between breathers to 
Anderson modes are available in |KA99-00t IAMM99j . 

In addition, the existence of breathers in inhomogeneous Klein-Gordon lattices 
(with disordered on-site potentials) has been proved by Sepulchre and MacKay 
|SM98l ISM97j for small coupling k. The proof is based on the continuation method 
previously introduced by MacKay and Aubry |MA94j for an homogeneous chain 
(method considerably generalized in [SM97j ). For k = the system reduces to 
an array of uncoupled non-identical anharmonic oscillators, and the simplest type 
of discrete breather consists of a single particle oscillating while the others are at 
rest. Under a nonresonance condition |SM98l [SM97j . this solution can be continued 
to small values of k (in most cases at fixed frequency) using the implicit function 
theorem, yielding a spatially localized solution. 

In this paper we provide complementary mathematical tools for studying time- 
periodic oscillations (not necessarily spatially localized) in inhomogeneous infinite 
lattices. The theory is developed in a very general framework, and applied to 
breather bifurcations in inhomogeneous Klein-Gordon lattices as lattice parameters 
and breather frequencies are varied. We start from a general Klein-Gordon lattice 
with Hamiltonian 

H= -Y^l + DnV{AnXn) + ^{Xn+l-Xnf (2) 

n=— oo 

(case p = of ([T]) with more general inhomogeneities). The potential V is assumed 
sufficiently smooth in a neighbourhood of with V^'(O) = 0, V"{0) = 1. The general 
theory is a priori valid for small inhomogeneities and small amplitude oscillations. 
In particular, in our application to system ([2]) we assume M„, D„, An, Kn to be 
close (uniformly in n) to positive constants. However, considering an example of 
Klein-Gordon lattice with a mass defect, we check using numerical computations 
that our tools remain applicable up to strongly nonlinear regimes, and sometimes 
for a large inhomogeneity. 

Our analysis is based on a centre manifold reduction and the concept of spatial 
dynamics. This concept was introduced by K. Kirchgassner |Kir82j for nonlinear 
elliptic PDE in infinite strips, considered as an (ill-posed) evolution problem in the 
unbounded space coordinate, and locally reduced to a finite-dimensional ODE on 
an invariant centre manifold. This idea was transposed to the context of travelling 
waves in homogeneous infinite oscillator chains by looss and Kirchgassner [IKOOj . 
and centre manifold reduction has been subsequently applied to the analysis of trav- 
elling waves and pulsating travelling waves in different one- dimensional homogeneous 
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lattices |looOOl [J505l llJOSl EmEI IPROSIUPOG] . Indeed, looking for travelling waves 
in an oscillator chain yields an advance-delay differential equation (a system of such 
equations in the case of pulsating travelling waves), which can be reformulated as an 
infinite-dimensional evolution problem in the moving frame coordinate, and locally 
reduced to a finite-dimensional ODE under appropriate spectral conditions. 

In reference |Jam01j . one of us has proved the existence of breathers in Fermi- 
Pasta- Ulam (FPU) lattices using a similar technique in a discrete context. The 
dynamical equations for time-periodic solutions were reformulated as an infinite- 
dimensional recurrence relation in a space of time-periodic functions, and then lo- 
cally reduced to a finite-dimensional mapping on a centre manifold, where breathers 
corresponded to homoclinic orbits to 0. A general centre manifold theorem for 
infinite-dimensional maps with unbounded linearized operator has been proved sub- 
sequently |Jam03j and has been used to analyze breather bifurcations in diatomic 
FPU lattices |JN04l IJKOTj and spin lattices |Nob04] . More generally, the dynam- 
ical equations of many one-dimensional lattices can be reformulated as infinite- 
dimensional maps in loop spaces as one looks for small amplitude time-periodic 
oscillations ( |Jam03j . section 6.1). 

As shown in the present paper, the centre manifold reduction theorem readily 
applies to homogeneous Klein-Gordon lattices, where M„ = m, Dn = d, An = a, 
Kn = /c in ([2]) and m,d,a,k > 0. This reduction result rigorously justifies (in 
the weakly nonlinear regime) a formal one-Fourier mode approximation previously 
introduced in reference |BCKRBW00j . The equations of motion read 

m—^ + daV'{aXn) = k{xn+i-2xn + Xn-i), n G Z. (3) 

Looking for time-periodic solutions (with frequency a;) and setting x„(t) = Xn{ujt), 
([3]) can be formulated as an (ill-posed) recurrence relation {xn+i,Xn) = F{xn,Xn-i) 
in a space of 27r-periodic functions. Using the theorem of reference | Jam03] . one 
can locally reduce the problem to a finite-dimensional mapping on a centre manifold 
whose dimension depends on the frequency u. More precisely, equation linearized 
at x„ = admits solutions in the form of linear waves (phonons) with Xn{t) = 
A cos {qn — Uqt) , whose frequency satisfies the dispersion relation 

mu;^ = a'^d + 2k{l- cosq). (4) 

The frequencies Uq lie in a band [uJmin,i^max] with utmin > 0. In the nonlinear 
case, the dimension of the centre manifold depends on how many multiples of u 
belong to (or are close to) the phonon band. When u ~ uJmax or ~ ujmin (with 
no additional resonance), the centre manifold is two-dimensional if solutions are 
searched even in time, which reduces ([3]) locally to a two-dimensional reversible 
mapping on the centre manifold. For appropriate parameter values, this map admits 
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small amplitude homoclinic solutions to corresponding to breather solutions of 
1^. Breather solutions in this system have been proved to exist by MacKay and 
Aubry |MA94] for small values of the coupling parameter k. Known regions of 
breather existence are considerably extended here, since we prove the existence of 
small amplitude breathers for arbitrary values of k in some cases and for frequencies 
close to the phonon band edges (see theorem [71 pJHTI). In addition we prove the 
existence of "dark breather" solutions, which converge towards a nonlinear standing 
wave a.s n ±oo and have a much smaller amplitude at the centre of the chain. 
These solutions correspond to heteroclinic orbits of the reduced two-dimensional 
map. 

Furthermore, we extend this analysis to the case when small lattice inhomo- 
geneities are present. The dynamical equations of the inhomogeneous system 
take the form 

(Px 

1^ ^ ^""^^ V(A„X„) = Kn (X„+1 - Xn) - Kn-l{Xn - Xn-l), U E Z, (5) 

and time-periodic solutions can be obtained as orbits of a nonautonomous map 
{xn+i,Xn) = F{Xn,Xn,Xn-i), whcrc the noucoustaut lattice parameters are embed- 
ded in a multicomponent parameter A„. Fixing M„ = m + m„, Dn = d + dn, 
An = a + an, Kn = k + kni we consider the case when constant lattice parameters 
'm,d,a,k > are perturbed by uniformly small sequences {mn)nGZ, {dn)n&z, {an)nez, 
{kn)n&- We prove (see theorem [HI pJT^ that small amplitude time-periodic solutions 
with frequencies close to oJmm or ^max are determined by a two-dimensional nonau- 
tonomous map. Moreover, we generalize this reduction result in the case when 
several multiples of u) are close to the band [uJmin^^max], which yields a higher- 
dimensional reduced problem (see theorem HI p ITOl) . 

In fact we prove this type of reduction result in a very general framework, for 
infinite-dimensional mappings with small nonautonomous perturbations, considered 
in a neighbourhood of a non-hyperbolic fixed point, or close to a bifurcation. The 
linear autonomous part of the map must satisfy a property of spectral separation 
(see theorem [T] pIHj), but a large number of one-dimensional lattices with finite- 
range coupling fall within this category. We obtain a direct proof of the reduction 
result by observing that any nonautonomous mapping Un+i = F{Xn, Un) can be seen 
as a projection of an extended autonomous mapping, to which the centre manifold 
theorem of reference |Jam03] can be applied under appropriate assumptions. The 
centre manifold of the extended map is infinite-dimensional, but this case is also 
covered in |Jam03j . The reduced nonautonomous mapping for the original system 
can be interpreted as a projection on a finite-dimensional subspace of the extended 
autonomous mapping restricted to the invariant centre manifold. 

We use this reduction result to analyse the case when equation ([5]) presents a 
mass defect at a single site, all other lattice parameters being independent of n. In 
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that case, the hnearized problem admits a spatially localized mode (usually denoted 
as an impurity mode), and a nonlinear continuation of this mode can be computed 
|FPM94j . corresponding to a Lyapunov family of periodic orbits. Klein-Gordon sys- 
tems with a coupling defect or a harmonic impurity in the on-site potential share 
similar characteristics |CPAR02] . as well as nonlinear lattices with a different type 
of nonlinearity |KZK97j . In addition to this simple localization phenomenon, sin- 
gle impurities can have more complex effects in a nonlinear system. Indeed it is a 
common feature to observe a complex sequence of tangent bifurcations between (de- 
formations of) site-centered and bond-centered breathers in some neighbourhood of 
the defect as the strength of an impurity is varied jCK04l [KZK97j . Using numerical 
computations we show some examples of such bifurcations in the present paper, as 
one varies the strength of a mass defect in system . From a physical point of view 
it is quite important to understand how a local change in the lattice parameters 
modifies the set of spatially localized solutions. For example, this could contribute 
to explain how a mutation at a specific location of an homogeneous sequence of 
(artificial) DNA would modify the structure of fluctuational openings |KRBCU04] . 

This paper provides a qualitative explanation of such tangent bifurcations, which 
reveals also very precise quantitatively when compared with numerical simulations 
of the Klein-Gordon model. According to the previously described reduction theo- 
rem, for a small mass defect of size e, small amplitude breather solutions of ([5]) with 
frequencies below (and close to) ujmin are described by a two-dimensional nonau- 
tonomous mapping Vn+i = f{vn,oj) + eg{n,Vn,oj,e). Here we only study the prin- 
cipal part of the reduced mapping as (t>„,u;,e) ~ (0,co'mm5 0). We show that this 
truncated reduced map admits an homoclinic orbit to (corresponding to an ap- 
proximate breather solution for the oscillator chain) if, for e = 0, the image of the 
unstable manifold of under a certain linear shear intersects its stable manifold. 
The linear shear is 0(e)-close to the identity. When the on-site potential is soft, 
these manifolds have very complicated windings characteristic of homoclinic chaos, 
hence the set of their intersections changes in a complex way as the linear shear 
varies, or equivalently as one varies the mass defect. This phenomenon explains the 
existence of the above mentioned tangent bifurcations, at least for small defect sizes, 
and for small amplitude breathers with frequencies close to the phonon band. In 
addition, we show (by comparison with direct numerical simulations of the Klein- 
Gordon model) that this picture remains valid quite far from the weakly nonlinear 
regime. Let us note that, to obtain an exact solution of ([5]) from an orbit of the 
truncated map, it would be necessary to control the effect of higher order terms 
(with respect to f„, u — Umin, e) present in the full reduced mapping and prove the 
persistence of this solution. This problem is not examined here from the analytical 
side, but we compare instead numerically computed solutions of (ED with approxi- 
mate solutions deduced from the truncated map. The very good agreement leads 
us to conjecture that most of the homoclinic bifurcations existing for the truncated 
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problem persist for the full reduced system. 

Lastly we consider the more general case when system admits a finite number 
of defects, i.e. perturbations m„, dn, an, kn have a compact support (as above these 
perturbations are assumed to be small, of order e). We show that the approach 
developped for a single impurity can be extended to this case (see lemma 171 p HSll . 
where the linear shear is replaced by a more general linear near-identity transforma- 
tion A^. The linear transformation provides a useful tool for studying breather 
bifurcations in Klein-Gordon lattices with a finite number of impurities, as for the 
single impurity case that we have analyzed in detail. By computing the principal 
part of A^ as e is small and frequencies are close to Umin, we show that the effect 
of the parameter sequence on the set of small amplitude breather solutions should 
mainly depend on weighted averages of the defects values. 

The outline of the paper is as follows. Section [2] presents the centre manifold 
reduction theory for time-periodic oscillations in weakly inhomogeneous nonlinear 
lattices. We treat the case of Klein-Gordon lattices in detail in sections 12.11 and 
12. 3[ and formulate the reduction theory in a much more general setting in section 
12.21 Section [S] concerns spatially homogeneous Klein-Gordon lattices. Existence 
theorems for small amplitude breather and dark breather solutions are deduced from 
the dynamics of two-dimensional reversible maps on invariant centre manifolds. The 
case of weakly inhomogeneous Klein-Gordon chains is considered in section HI where 
the truncated reduced map is analyzed for a finite number of defects. A geometrical 
condition for the existence of homoclinic orbits to is derived in section 14. 2^ and 
some homoclinic bifurcations are studied in detail in section 14.11 for a single mass 
defect. In the latter case, breather solutions are numerically computed in section [5] 
and the results are successfully compared with our analytical findings. 

2 Reduction result for small inhomogeneities 

In this section we consider system ([5]) in the limit of small inhomogeneities. We 
show that all small amplitude time-periodic solutions are determined by a finite- 
dimensional nonautonomous map, whose dimension depends on the frequency do- 
main under consideration. For this purpose we reformulate ([S]) as a map in a loop 
space, perturbed by a small nonautonomous term (section 12.11) . Then we prove 
in section 12.21 a general centre manifold reduction theorem for infinite-dimensional 
maps with small nonautonomous perturbations. This result is based on the centre 
manifold theorem proved in reference |Jam03j for autonomous systems. Our general 
result is applied to the inhomogeneous Klein-Gordon lattice, which yields the above 
mentioned reduction result (section [23]) • 
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2.1 The Klein-Gordon system as a map in a loop space 

We set Xn{t) = yniuj{k/my^'^t) in equation ([5]), where ?/„ is 27r-periodic in t (hence 
Xn is time-periodic with frequency u{k/mY/'^). The constant a > being fixed, we 
also define V{x) = a~'^V{ax). Equation ((51) becomes 

(6) 

where = a'^d/k and l + e„ = (1 + + 1 + = (1 + ^)(1 + ^)/(l + 

7„ = ^, 1 + = (1 + ^)/(l + The sequences (e„)„6z, (^?n)nez, (7n)nez 
('^n)nez will be assumed sufficiently small in iao{Z), where ^oo(^) is the classical 
Banach space of bounded sequences on Z, equiped with the supremum norm. To 
simplify the notations, we shall drop the tilde in the sequel when referring to the 
renormalized potential V. Moreover we shall use the shorter notations {e} when 
referring to sequences (e„)„gz- 

To analyze system ([6]) we use the same approach as in reference jJam03j for spa- 
tially homogeneous systems. We reformulate ([6]) as a (nonautonomous) recurrence 
relation in a space of 27r-periodic functions of t, and locally reduce the (spatial) 
dynamics to one on a finite-dimensional centre manifold. We restrict our attention 
to the case when i/n is even in t in order to deal with lower-dimensional prob- 
lems. More precisely, we assume ?/„ G for all n G Z, where H!^ = { ?/ G 
ifpg^(0, 27r), y is even} and ifpg^(0,27r) denotes the classical Sobolev space of 2it- 
periodic functions (if°g^(0, 27r) = L^g^(0, 27r)). 

Since our analysis concerns small amplitude solutions and small inhomogeneities, 
the first step consists in studying the linearized system at = when e„, 7^, k„ 
are fixed equal to 0. In that case equation yields 

uj'^^^^ + ^'^yn = yn+i-'2yn + yn-i, U E Z. (7) 

Now we rewrite the problem as an infinite-dimensional linear mapping. For this 
purpose we introduce F„ = (?/„_i,?/„) G D, where D = x H"^. Equation can 
be written 

= Yn, n G Z, (8) 
where ^ 

A^{z,y) = (^y , uj^'^ + + 2)y - z) (9) 

and equation ([H]) holds in X = H"^ x H^. The operator : D d X X is 
unbounded in X (of domain D) and closed (we omit the additional parameter f2 in 
the notation A^). 
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The spectrum of consists in essential spectrum at the origin and an infi- 
nite number of eigenvalues Up, {p > 0) depending on a;, and satisfying the 
dispersion relation 

+ (cjV - _ 2)(j + 1 = (10) 

(it follows that is either real or has modulus one). Equation ( ITOl) is directly 
obtained by setting Hn = cr"' cos (pt) in equation ([7]). The invariance a in 
f llOl) originates from the invariance n —>■ —n in ((Tj). In the sequel we shall note Up 
the solution of ffTUl) satisfying |crp| > 1 and Imo"p < 0. Clearly ap is real negative 
for p large enough and limp^+ooO'p = — oo. Moreover accumulates at cr = as 
p — s> +00. It follows that the number of eigenvalues of A^^ on the unit circle is finite 
for any value of the parameters uj,Q. 

In addition, the eigenvalues cXp, defined by ( [TUl) lie on the unit circle when 
n < up < (4 + f2^)^/^. This property as a simple interpretation. Multiplying fIlOp by 
a~^, setting a = e*"^ and cUg = u p (k/mY^'^, one finds the usual dispersion relation 
(jl]). Consequently, if ujp{k/mY^'^ lies inside the phonon band [uJmin,^max] for some 
p G N, then A^^ admits a pair of eigenvalues e^""^ on the unit circle determined 
by the dispersion relation (jl]). This condition on uj is equivalent to prescribing 
n<Lup< (4 + fi2)i/2_ 

Now let us describe the spectrum of A^^ near the unit circle when > is 
fixed and u is varied. As we shall see, the number of eigenvalues of A^^ on the 
unit circle changes as uj crosses an infinite sequence of decreasing critical values 
cji > ci;2 > • • • > 0. Small amplitude solutions of the nonlinear system bifurcating 
from yn = will be found near these critical frequencies. 

We begin by studying the evolution of each pair of eigenvalues Cp, ap^ as u 
varies. Firstly, one can easily check that uo, a^^ are independent of to, real positive 
and lie strictly off the unit circle. 

Secondly we consider the case p > 1. For u > y/4+Tfi/p, (Jp^a^^ are real 
negative and lie strictly off the unit circle. When uj decreases, they approach the 
unit circle and one has Up = = —1 for cu = V4 + fi^/p (this corresponds to 
a frequency Uq at the top of the phonon band, for a wavenumber q = n). At this 
critical parameter value, cXp = —1 is a double non semi-simple eigenvalue of A^. 
For Q/p < UJ < y/4+Ti^/p, ap,ap^ lie on the unit circle, and approach +1 as u 
decreases. One has CTp = a'^ = 1 for = Q/p, and then +1 is a double non semi- 
simple eigenvalue of A^ (this corresponds to a frequency Ug at the bottom of the 
phonon band, for a wavenumber q = 0). For lu < fl/p, ap, ap^ are real positive and 
lie strictly off the unit circle. 

Now let us qualitatively describe the evolution of the whole spectrum of A^. 
When UJ > v4+T2^ the spectrum of lies strictly off the unit circle (both inside 
and outside). When uj decreases, the eigenvalues dp approach the unit circle for 
all p > 1. As the first critical value uji = a/4 + is reached, the eigenvalues 
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ai,crf^ collide and yield a double (non semi-simple) eigenvalue ai = —1, while the 
remaining part of the spectrum is hyperbolic. When 00 is further decreased, two 
different situations occur depending on the value of fl. 

For Q > 2/V3, 0"! , cr^ ^ are the only eigenvalues on the unit circle for n < cu < 
y/^+Tfi. One has cxi = ajf ^ = 1 at the second critical value uj2 = fl. When to is 
further decreased, (Ti,(7f^ are real positive and lie strictly off the unit circle. One 
has (T2 = (72^ = — 1 at the third critical value uj^ = v4 + ^/2 < fl. The situation 
is sketched in figure [H 




Case uj > ui Case L02 < ^ < ^^i Case 003 < uj < UJ2 

Figure 1: Spectrum of near the unit circle as u is varied, in the case 
The unbounded part of the spectrum on the negative real axis is not shown. The 
arrows indicate how the eigenvalues have moved from their positions in the previous 
graph, after lj has been decreased. 

The case < 2/^3 is different, since ai , ^ are the only eigenvalues on the unit 
circle in the smaller frequency range v4+T2^/2 < u < Indeed one has 

(T2 = cr^^ = —1 at second critical value 102 = V4 + f2^/2 > il. For lu < V4 + i7^/2 
and u ~ a/4+TF/2 the spectrum of on the unit circle consists in two pairs 
of simple eigenvalues (Xi, ^, (72, o"^^. In the interval < to < + other 
eigenvalues may collide at —1 depending on the value of Q. The situation is sketched 
in figure [21 

In what follows we restrict our attention to the neighbourhood of critical fre- 
quencies uj ^ UJ2 with r2 > 2/-\/3, and ~ tui. This leads us to consider the small 
parameter /i defined by = uf + fi. As u equals one of the critical frequencies 
uJi,W2, the spectrum of A^ on the unit circle only consists in a double eigenvalue 
— 1 or +1. For u) ^ Ui or u)2, the above spectral analysis shows that the fixed point 
y = of ([H]) is hyperbolic when uj > uJi or uj < 002- In this case, when nonlinear 
effects will be taken into account, we shall see that the stable and unstable manifolds 
iy*(0), Vl^"(0) may intersect depending on the local properties of the anharmonic 
potential V, leading to the existence of homoclinic orbits to y = 0. 

System ([6]) will be analyzed in the limit of small amplitude solutions and for 
small parameters /i, {e}, {77}, {7}, {k}. The parameter space will be denoted as 
E = Rx (£00 (Z)) - AH parameters are embedded in a multicomponent parameter 
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Case Lo > LUi Case u!2 < ^ < Case uj^ < u < uj2 

Figure 2: Spectrum of A^^ near the unit circle as lu is varied, in the case n < 2/^3. 
The unbounded part of the spectrum on the negative real axis is not shown. The 
arrows indicate how the eigenvalues have moved from their positions in the previous 
graph, after u has been decreased. 

{A} = (fi, {e}, {rj}, {7}, {k.}) G E. In addition we denote by r„ the index shift in 
£oo(Z), i.e. {Tn{e}}k = en+k- 

Equation ^ can be rewritten in the form of a nonautonomous mapping in a 
function space. More precisely we have 

F„+i = LF„ + iV(y„,A„), nGZ, (11) 

where F„ = (?/n-i,2/n) = {zn^Vn) e D, Xn = {fi,en,r]n,'ln,fin) G L = A^^ (for 
i = 1 or 2) and N{z, y, A„) = ( , N2{z, y, A„) ), 

N2{z,y,\n) = (^fe„+/i(l+e„)) -j^+n^[{l+nr,){l+^n)-l]y+^^n{y-z)+W{y,r]r^,^n), 

Wiy, r^,^) = n'il + r]){V'[il + ^)y]-il + ^)y). (12) 

Equation (fTTl) holds in the Hilbert space X. The potential V is assumed sufficiently 
smooth {C'^~^^, with p > 5) in a neighbourhood of 0. It follows that : DxMP —>■ X 
is {k = p — 2 > 3) in a. neighbourhood of (Y, A) = 0. The operator consists in 
higher order terms as {Y, A„) ^ 0, i.e. we have A^(0, A) = 0, DyN{0, 0) = 0. 

We note that ( JTTl) is invariant under the symmetry TY = Y{- + 7i). Moreover the 
usual invariance under index shifts {Y} ti{Y} is broken by the inhomogeneity 
of the lattice, and replaced by the invariance {{Y}, {A}) — > {ti{Y} , ti{X}) . 

In the next section we prove a general centre manifold reduction theorem for 
maps having the form fill I) , under appropriate spectral conditions on L and for 
small nonautonomous perturbations {A} G E. This analysis relies on the reduction 
results proved in reference |Jam03j for autonomous maps. To simplify the proof, 
problem (fTTl) will be considered as a projection of a suitable autonomous mapping 
to which the centre manifold theorem can be directly applied. 
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2.2 Centre manifold reduction for nonautonomous pertur- 
bations of infinite-dimensional maps 

In this section we reformulate the situation of section 12.11 in a general framework, 
and prove a local centre manifold reduction result for problems of this type. This 
level of generality is relevant for nonlinear lattices, because the dynamical equations 
of many one- dimensional lattices can be reformulated as infinite-dimensional maps 
in loop spaces as one looks for small amplitude time-periodic oscillations. Indeed, if 
the coupling between sites has finite range (i.e. Xn is coupled to Xk for \n — k\ < p), 
then in general Xn+p can be obtained locally as a function of Xn+p-i, • • • , Xn-p using 
the implicit function theorem (for some examples see e.g. j Jam03] . section 6.1, or 
jJKOT] ). 

To work in a general setting, let us consider a Hilbert space X and a closed linear 
operator L : D G X ^ X of domain D, L being in general unbounded. We equip 
D with the scalar product {u, v)j^ = {Lu, Lv) ^ + (m, v)-^, hence D is a Hilbert space 
continuously embedded in X. 

We denote by W x V a neighbourhood of in D x M''' and consider a nonlinear 
map G C"=(WxV,X) {k > 2), where A^(Y', A) satisfies A^(0, A) = 0, DyN{0,0) = 0. 
We look for sequences (Yn)nez in ^ satisfying 

WneZ, Yn+i = LYn + N{Yn,Xn) inX, (13) 

where {A} = {\n)n& is a bounded sequence in V treated as a parameter. In what 
follows we shall note E = £oo(^, I^^) the Banach space in which {A} lies. Notice 
that y = is a fixed point of (|T3l) . 

We assume that L has the property of spectral separation, i.e. L satisfies the 
assumption (H) described below (in what follows we note o"(T) the spectrum of a 
linear operator T). 

Assumption (H): The operator L has nonempty hyperbolic (\z\ 1) and central 
(\z\ = 1) spectral parts. Moreover, there exists an annulus A = {z E C , r < \z\ < 
R} (r < 1 < R) such that the only part of the spectrum of L in A lies on the unit 
circle. 

The situation corresponding to assumption (H) is sketched in figure El Under 
assumption (H), the hyperbolic part ah of cr{L) is isolated from its central part o"c. 
In particular this allows one to split X into two subspaces Xc, Xh invariant under 
L, corresponding to cJc, an respectively. More precisely, Lh = L\x^ and Lc = L\x^ 
satisfy a{Lh) = ah and a{Lc) = ac- 

The invariant subspace Xc is called centre subspace, and Xh is the hyperbolic 
subspace. The subspace Xc is finite-dimensional when the spectrum of L on the 
unit circle consists in a finite number of eigenvalues with finite multiplicities (we do 
not need this assumption for the reduction theorem constructed here). 
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The spectral projection tt^ on the centre subspace can be defined in the following 
way (see e.g. |Kat66j ) 

T^c = 7^ [ {zI-L)-Uz-^[ {zI-L)-Uz, 

where C{r) denotes the circle of centre z = and radius r (see figure E])- One has 
TTc G C{X,D), Xc = C D and = LtIc, where C{X,D) denotes the set of 
bounded operators from X into D. In the sequel we note tt/i = / — tTc and Dh = i^hD- 
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Figure 3: Spectrum of L (dots), unit circle (dashed) and oriented circles C(r), C{R). 

Remark 1 Let us consider the situation of section \2.1\ and the linear operator L of 
equation 07]) . In the case uj = uj2 and Vt > 2/\/3, the spectrum of L on the unit 
circle consists in a double non semi-simple eigenvalue +1. Moreover, for uj = ui 
the spectrum of L on the unit circle consists in a double non semi-simple eigenvalue 
— 1. In both cases the associated invariant subspace X^ is spanned by = (cost, 0), 
Vy = (O,cost) and we have tCcY = ^{J^^ Y{t) cost dt) cost. In addition the unit 
circle is isolated from the remainder of the spectrum, since the latter is discrete and 
only accumulates at the origin and at — oo on the real axis. It follows that L satisfies 
assumption (H). 

Now we state the centre manifold reduction theorem in the general case. In the 
sequel we note = tTc Y, Y^ = Tih Y. 

Theorem 1 Assume that L has the property of spectral separation, i.e. satisfies 
assumption (H). There exists a neighbourhood QxAofOinDxE and a map 
^ e C'iX^ X A, Dh) (with ^(0, {A}) = 0, L'yc^/>(0, 0) = o) such that for all {A} G A 
the following holds. 

i) If {Y} is a solution of ( fl^) such that y„ G i7 for alln G Z, then Y^ = ipiY^, Tn{A}) 
for all n eT: and Y^ satisfies the nonautonomous recurrence relation in X^ 

VnGZ, y„V = /n(^n,{A}), (14) 



14 



where fn G C^{{Xc HQ) x A,Xc) is defined by 

/„(., {A}) = vr, (L + N{., A„)) o (/ + r„{A})). 

ii) Conversely, if {Y'^} is a solution of [14^ such that E Q for all n E Z, then 
Yn = Y^ + ,p{Y:i, T„{A}) satisfies 

Hi) IfL+N{., A) commutes with a linear is ometry T G C{X)r\jC{D) thenTip{., {A}) = 
^|J{., {A}) o T and Tf^., {A}) = /„(., {A}) o T. 

Properties i), ii) reduce the local study of ffT^ to that of the nonautonomous 
recurrence relation ( |T4l) in the subspace Xc. Note that the dependency of and the 
reduced map /„ with respect to sequences {A} is nonlocal. 

In what follows we give a simple proof of theorem [T] which relies on the fact 
that the nonautonomous mapping (fT3i) can be seen as a projection of an extended 
autonomous mapping, to which the centre manifold theorem proved in |Jam03] can 
be applied. This procedure will explain why the result of theorem [1] can be seen 
as a centre manifold reduction, since the reduction function ip will appear as one 
component of the function having the centre manifold as its graph. The reduced 
nonautonomous mapping (iMl) will be interpreted as a projection of the extended 
autonomous mapping restricted to the invariant centre manifold. 

Theorem [T] has been proved in reference [Jam03j in the case of an autonomous 
mapping, when the sequence {A} is absent or replaced by a simple parameter A G MP. 

To recover this autonomous case we introduce the additional variable Sn = 
Tn{A} G E. Note that for any fixed n ^ Sn denotes a bounded sequence in 
MP (to simplify the notations we use the symbol Sn instead of {Sn})- Given a 
sequence {A} E E we also note So{X} = Aq. Equation ffT^ can be rewritten 

Yn+l = LYn + N{Yn,6oSn), Sn+1 = T.Sn, (15) 

which consists in an autonomous mapping in X x E. 

In what follows we apply the theory of reference |Jam03] to system (1151) . As we 
shall see the corresponding centre manifold will be infinite-dimensional due to the 
second component of (fTSjl . The case of infinite-dimensional centre manifolds has 
been treated in reference |Jam03j . with the counterpart that theory is restricted to 
maps in Hilbert spaces. Consequently the first step is to search for Sn in a suitable 
Hilbert space instead of the Banach space E. For this purpose we consider the space 
of sequences 

h^, = {{u}/ukEC^, ||M||_i<+oo}, 

where = J2kez('^ ~^ ^'^)^^ W'^kW'^ ■ '^^^ space defines a Hilbert space 

equiped with the scalar product {{u},{v}) = Sfe^z + ^^)~^ ""fc ' ^fc' where ■ de- 
notes the usual scalar product on and || || the associated norm. For all n E Z we 



15 



now search for Sn in the space H = h^i fl (R^)^ consisting of real sequences in h_i. 
Note that E (Z H, the embedding being continuous. 

Since sequences in H may be unbounded and N(Y, .) is defined on a neighbour- 
hood V of A = in M^, we replace f lTSj) by a locally equivalent problem 

(F„+i,^„+i)=F(y„,^„) (16) 

where 

F{Y, S)={LY + iV(r, 7(5o5) ),nS), 

7 : IRP — > V is a C°° cut-off function satisfying ||7(x)|| < ||a;||, 7(0;) = x for ||a;|| < r, 
7(0;) = for ||a;|| > 2r, r being chosen small enough (with 5(0, 2r) C V). 

Problem (fTBj) consists in an autonomous mapping in X x H. In order to apply 
the centre manifold theorem of reference [Jam03] we need to study the spectrum 
of DF{0) = L X n. One has clearly a{DF{0)) = a{L) U o-(ri), where a(ri) is 
determined in the following lemma. 

Lemma 1 The spectrum cr(ri) of Ti : H H consists of the unit circle. 

Proof. Consider the complexification h-i of H. Given a sequence {/} G and 
2; G C, we look for {u} G /i-i satisfying 

{zl-n){u} = {f}. (17) 

Equation ( ITTl) can be solved in a simple manner using Fourier series. Recall that 
the periodic Sobolev space ifpg^(0,27r) can be defined as the set of functions in 
L^(M/27rZ, C^) whose Fourier coefficients form a sequence in hi, where 

hi = { {u} /ukE CP, J] + ll^'^ll^ < +^ ^• 

fcez 

In the same way its dual space H~Jj,{0,27i) is isomorphic to h_i, where the iso- 
morphism C : H~J^{0,2n) h_i is again given by CniT) = ^(T, e~*"*) for all 
T G H~J^{0,27i). In addition one has the useful property riC(T) = C(e~**T) for 
all T G ifj^J.(0, 27r). Now return to equation (ITTI) and consider T = C^^dw}) and 
S' = C^^{{f}). One obtains the equivalent problem in H~J^{0,2Tr) 

{z-e-'')T = S. (18) 

If \z\ 7^ 1 then (|TS|) has the unique solution T = {z — e~**)~^S', hence z ^ o"(ti)- If 
2; = e*^ is chosen on the unit circle, T = 27r5_e is a solution for S" = 0, corresponding 
to an eigenvector {u} = {e*"^} of ri. 

□ 
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As it follows from lemma [Tj cr(i5-F(0)) consists of the union of cr(L) with the 
unit circle. Consequently DF{0) has the property of spectral separation, i.e. the 
hyperbolic part of its spectrum is isolated from the unit circle. Moreover the centre 
subspace of DF{0) is simply Xc x H. With these spectral properties at hand, we now 
apply the centre manifold theorem of reference |Jam03] which states the following. 

Theorem 2 There exists a neighbourhood Q x A of (Y, S) = in D x H and a map 
^ G C^(Xc X H, Dh) (with ^{Q, 0) = 0, ^^^(0, 0) = such that the manifold 

M = {{Y,S) e D X H /Y = Y^ + S),Y^ eX^} 

has the following properties. 

i) M is locally invariant under F, i.e. if ("K, S) E M f^ {Vl x K) then F{Y, S) E M. 

a) If {{Y, S)} is a solution of ( fl^) such that {Yn, Sn) & Q x A for all n G Z, 
then (Yn, Sn) G M. for all n E Ij (i.e. Y^ = ip{Y^,Sn)) and {Y^,Sn) satisfies the 
recurrence relation in Xc x H 

Y^_^_^ = f{Y^, Sn), Sn+l = TiSn, (19) 

where 

f{Y% S) = TT, [L + iV(., i{5,S))]{Y^ + S)). 

Hi) Conversely, given a solution {{Y'', S)} of (13) such that (K^, Sn) G x A for all 
n G Z, consider Yn = Y^ + ip(Y^, Sn). Then (Yn, Sn) defines a solution of ( fl^) lying 
on M. 

iv) IfL+N{., A) commutes with a linear is ometry T G C{X)r\C{D) thenTiplY'^, S) = 
iIj{tY\ S) and Tf{Y^, S) = f(TY^, S). 

The manifold M. is called a local centre manifold for ffTBj) . It is locally 
invariant under F (as stated by property i) ) and the linear isometrics of ffTUl) . 
Property iv) expresses the invariance of TVI under the linear isometry T x I of ( fT6l) . 

Now the proof of theorem [T] follows directly from theorem [21 Since E is contin- 
uously embedded in H, defines a C'^ map from X^x E into D^. In theorem [T] 
we choose A as a ball of centre Q m. E such that A C A and 7 = / on A. Then 
problems ( fT3l) and ( |T6i) are equivalent for all {A} G A, with Sn = '^^{A}, and proper- 
ties i)-ii)-iii) of theorem [1] are directly deduced from properties ii)-iii)-iv) of theorem 
[21 In addition, since (0, t„{A}) is a solution of (fTUI) for all {A} G A it follows 
'?/'(0, r„{A}) = (by property ii) of theorem [2]), and consequently ip(Q, {A}) = 0. 
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2.3 Application to the Klein-Gordon lattice 
2.3.1 Reduction result 

In this section we apply the reduction theorem[T]to the inhomogeneous Klein-Gordon 
lattice considered in section 12.11 We recall that the inhomogeneous system ([6]) has 
been reformulated as a nonautonomous map in a loop space given by expression fill I) . 
All parameters (sequences of heterogeneities and frequency shift /x) are embedded in 
the multicomponent parameter {A} = (/i, {e}, {r/}, {7}, {k}) E E = M. x (£oo(^))^- 
The problem has exactly the general form (JT3l) (in a particular case when the first 
component of {A} is constant) and consequently the reduction theorem [1] can be 
applied to ffTTl) . This yields the following result for the original system ([U]). 



Theorem 3 Fix uJ^ = uJl + ^ in equation where uo^ = v4+TF or Uc = (in 
that case we further assume > 2/^/3). There exist neighbourhoods U, V and W of 
in H^, E and M respectively, and a C'^ map : M?xE (with 0(0, {A}) = 0, 

D(f){0, 0) = 0) such that the following holds for all {A} G V. 

i) All solutions of ^ such that Un &U for all n E'L have the form 

yn{t) = PnCOSt + Hn{t), 

where Hn = 0(/3n-i, /5n, t„{A}). For lOc = ^, Pn satisfies a recurrence relation 

Pn+l - 2Pn + Pn-1 = RniPn-1, Pn, {A}) (20) 

where Rn '■ x V ^ is . The principal part of Rn reads 

Rn{a,P,{X}) = {n'r]nil + ^n)-{n' + fj.)en + n'^n-lJ^)P + ^^n{P-a)+BP^ + h.o.t., 

(21) 

B = ^{V^'\0)-l{V^'\0)r). (22) 



For ujc = a/4 + one has 

Pn+l + 2/?„ + Pn-1 = RniPn-1, Pn, {A}), (23) 

with 

Rn{a,P,{\}) = (fiV(l+7n)-(4 + fi'+/x)e„ + fiV-/^)/?+'^n(/?-«) + S/3=^ + /i.o.i, 

(24) 

B = ^{ V^'\0) + iV^'Wij^^, - 2) ). (25) 

In both cases, higher order terms in Rn are 0(||(q;,/3)||'^ ||{A}||£; + ||(a;,/3)||^) and 
non-local in {A}. 
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ii) If Pn is a solution of problem /[2U\) or (d^j (respectively for Uc = ^ or Uc = 
\/4 + il'^), such that Pn £ W for alln G then yn{t) = /3„ cost + r„{A}) 
satisfies equation 

Hi) The functions and Rn have the following symmetries 

-P, {A}) = T0(a, p, {A}), i?„(-«, -P, {A}) = -Rn{a, P, {A}), 
where T denotes the half period time shift [T(f){.)]{t) = [</>(.)] (t + vr). 

It is straightforward to check that system has the reduction properties i) and 
ii) described above since the equivalent system (fTTl) satisfies properties i) and ii) of 
theorem [1] (see remark [Tlp JT^ . However there remains to compute the exphcit forms 
(j2Ti) and (IMl) of the recurrence relations. These expressions do not simply corre- 
spond to the two-dimensional mapping ffT^ rewritten as a second order recurrence 
relation. In addition we rewrite ffT^ in normal form, i.e. we perform a polynomial 
change of variables which simplifies (fT^ by keeping only its essential terms. These 
computations will be the object of the next three sections. Property iii) is equivalent 
to property iii) of theorem [U where the symmetry T is the half period time shift 
which satisfies T\x^ = —I- 

A possible way of computing the reduced recurrence relation would be to consider 
the equivalent autonomous mapping ffT^ and use a classical computation scheme 
for centre manifolds of autonomous systems (see e.g. |Van89] for a description of the 
method). The first step consists in computing the Taylor expansion of the reduction 
function up to a given order. This can be done using a nonlocal equation for 
(obtained by expressing F„ in ffT^ as a function of N{Yn, SqSu)) and computing the 
Taylor coefficients of ip by induction (see |Van89j ). The second step is to compute 
the reduced recurrence relation ( |T9i) which is completely determined by if). 

In the next three sections we shall use a different method yielding simpler compu- 
tations. Firstly we compute the recurrence relations fl2T|) and (12^ in the autonomous 
case {e} = {rj} = {7} = {n} = 0, using the method of reference | Jam03j . Then, 
using a symmetry argument, we deduce how the leading order part of the reduced 
equation is modified by the nonautonomous terms of (fTTl) . 

To end this section we point out a generalization of theorem [SI As it follows 
from the analysis of section 12.11 the dimension of the centre space X^. of is 
twice the number of multiples of cu lying within the band [fi, (4 + 1]^)^/^]. More 
precisely, if < tup < (4 + fi^)^/^ for p E {po, . . . ,pi}, with no additional multiples 
entering the band, then the centre space is spanned by the corresponding Fourier 
modes (cos(pt), 0), (0, cos(pt)). As above the following reduction result follows from 
theorem [H 

Theorem 4 Consider cUc > such that tOcP G [fi, (4 + fi^)-*^/^] for all integers p G 
{po, . . . ,pi} , with no additional multiples in this interval. Fix lo^ = ujI + ^ in 
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equation ^ and note N = pi — pq + 1. Consider the subspace He of spanned 
by the N Fourier modes cos{pot), . . . , cos(pit) and its complementary subspace 
consisting of orthogonal Fourier modes. There exist neighbourhoods U, V of in 
H"^, E respectively, and a map : x E (with 0(0, {A}) = 0, 

-D0(O, 0) = Oj, such that for all {A} G V, all solutions of ^ such that yn EU for 
all n ^ Z have the form 

pi 

yn{t) = [/^i'^ COS ] + '^(/^i-i' /^^^ • • • ' f^n-i r.M)- (26) 

P=PO 

Moreover, all small amplitude solutions of ^ are determined by a finite- dimensional 
recurrence relation obtained by projecting ^ on He and using the ansatz ^2E)- 

In the following sections 12.3.21 12.3.31 and 12.3.41 we compute the explicit forms of 
the reduced recurrence relations given in theorem [31 

2.3.2 Homogeneous case near the lower phonon band edge 

In this section we restrict our attention to the case when > 2/^/3 and to ^ UJ2 = ^■ 
We consider the autonomous case when {e} = {rj} = {7} = {k} = 0. Equation (ITTl) 
now reads 

= Lr„ + iV(F„,/x), nez (27) 

where L = is given by ([9]) and 

72 

with W{y) = ^^iV'{y) - y). System ([27]) is a reformulation of the equations of 
motion for the homogenous Klein-Gordon lattice 

u^-^ + Q^ V'iyn) = yn+i - 2yn + yn-i, neZ. (28) 

As in the nonautonomous case f[TT|) . system ([TTj) is invariant under the symmetry 
TY = Y{- + n). Moreover, the invariance y„ y^n of ([2H]) implies that ([TT]) is 
reversible with respect to the symmetry R{z, y) = {y, z), i.e. if Yn is a solution then 
also RY-n- In other words, if Y and [L + N{.,fi)]{RY) are in some neighbourhood 
of in D one has {L + A^(., yu) o R^Y = Y. Lastly, due to the existence of the addi- 
tional symmetry T, it is worthwhile to notice that TR defines an other reversibility 
symmetry. 

In what follows we use the notations introduced in section 12. 2[ We recall that 
the spectrum of L = Aq on the unit circle consists in a double non semi-simple 
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eigenvalue +1, and the associated two-dimensional invariant subspace X^. is spanned 
by the vectors K = (cost,0), Vy = (O,cost), with 



1 
-1 2 



in the basis {Vz, Vy). For fi in some neighbourhood A of 0, (|2711 admits a C'^ two- 
dimensional local centre manifold M.^ <Z D (which can be written as a graph over 
Xc), locally invariant under L + N{.,fi) (see |Jam03j . theorem 1 p. 32). One can 
write 

= {Y e D / Y = aVz + bVy + tPia, b, /x), (a, b) E }, (29) 

where ^ E C"=(M2 x A,Dh) and tp{a,b,fi) = 0{\\{a,b)f + \\{a,b)\\\fx\). Moreover, 
Ai^ is invariant under T and R (see |Jam03j . theorem 2 p. 34 and section 5.2). 

In the sequel we use the notations P*{y) = ^ J^^ y{t) cost (it, PcV = P*{y) cost 
and H"^ = {y E / P*{y) = }. The spectral projection tTc on reads Trdz, y) = 
{PcZ, Pcy) and we have Dh = Hi x Hi. 

Since A^^ is invariant under R and Vz, Vy are exchanged by R, we have the 
symmetry property Rip{a,b,fi) = ip{b,a,fi). Consequently, the function has the 
form 

ip{a, b, fi) = {(p{b, a, fi), (p{a, b, fi)) (30) 

with ip E C^(R'^ X A, Hi). Since A^^ is invariant under T and Tjx, = — / we have 
in addition 

T(p{a,b,iJ,) = (p{-a,-b,iJ,). (31) 

For yU ^ 0, the centre manifold contains all solutions Yn of ( |271) staying in 
a sufficiently small neighbourhood of F = in D for all n G Z. Their coordinates 
(a„, bn) on A^^ are thus given by a two-dimensional mapping which determines all 
small amplitude solutions when /i ^ 0. The reduced mapping is given by 

where 

( 6 ) " ( -a + 26 + r(a,6,/i) ) ' ^^^^ 
r{a,b,iJ,) = -iJ,b + P*W{bcost + ip{a,b,i2)). (34) 
One obtains equation flH^ using the fact that 

2:„ = a„ cost + v?(fen,an,Ai), 2/n = &n cost + V9(a„, 6„,/i) 

for Yn = {zn,yn) ^ Ai/i and applying P* to equation 1^7} (one has P*v9 = and 
P*o^ = -P* on Hi). 
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Since the reduced mapping inherits the symmetries of (1271) | Jam03j . commutes 
with T\x^ = —I and thus r(— a, — &, /i) = —r{a,b,fi). Moreover, is reversible with 
respect to the symmetry R{a, b) = (b, a), i.e. (/^ o = I. This yields the identity 

r(a, b, /i) = r(— a + 2b + r{a, b, /i), b, fi). 

This imposes the following structure for the Taylor expansion of r at (a, b, fi) = 

r(a, 6, /i) = + cifo^ + C2ab^ - ^C2a% + 0(|6| (|a| + + |6| (|a| + Vl), (35) 

where coefficients ci, C2 have to be determined. Note that r(a, 0, /x) = (see |Jam03j 
p. 53 for details). 

For determining the unknown coefficients of fl35p . we first compute the leading 
order terms in the Taylor expansion of ip at {a,b,^) = 0. This can be done using 
the fact that is locally invariant under L + N{.,fi) (see |Jam03j . theorem 1 p. 
32). For (a, 6) f« 0, this yields 

7ih[L + N{.,iJ,)]{{a,b) cost + ijj{a,b,n)) = ^^(/^(a, 6), /i) (36) 

or equivalently 

ip{-a + 2b + r{a, b, /i), 6, /i) = ip{a, 6, /x), (37) 
(p{b, —a + 26 + r(a, b, fi), fi) = 

{n^^ + 2 + n^)^{a, b, fi) - ^{b, a, fi) + {l- Pc) [fi^ + W] {b cos t + (^(a, b,fi)). 

(38) 

Thanks to the symmetry property fl37|) . the Taylor expansion of at order 2 takes 
the form ^ 

V?(a, 6, yu) = ifoiibfi - -ipii^a^ + v^iioo^ + V5o2o&^ + h-O-t- (39) 

By an identification procedure we now compute the coefficients ippqr in ( l39l) . using 
( l38i) and the expansion 

W{y) = n\ V(3)(0) + ^\/(^)(0) + 0{y') ). (40) 
Identification at order bfi gives 

(^ + l)v^on = 0, 

hence ipou = since v^oii ^ -f^^- Identification at order ab leads to 

1 £ 

^02o = -^{^^-j^ + 2 + n')^uo (41) 
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and identification at order 6^ yields 

72 



¥^110 + iS^^^ - 2 + n^)V020 = -^fiV(3)(0) cos^t. (42) 



Reporting fBTj) in fj^ gives 

,cP ... 2 



and consequently 



^110 = ^^^'^0) (1 + \ cos (2t)), ^020 = ^^(=^^0) (-^ - ^ + - ^) (2^))- 



As a conclusion, we obtain 

(43) 



¥p(a,6,/x) = ^\/(3)^o)(l + |cos(2t))(a6-ia2) 



+ i\/(=')(0) (-1 - J, + (1 - ^) cos (2t)) 62 + h.o.t. 

We now compute the two-dimensional mapping giving the coordinates (a„, h„ 
of the solutions on M.^. Equation flH^ can be written 



a„+i = 6„, - 26„ + 6„_i = r(6„_i, 6„, /x). (44) 

Using (IMl) . (HHjl and (jlD]) yields in equation 



ci = ^fiV(^)(o) - (^ + ^^^^)^(\/(^^o))^ c, = f^iy^'\Q)?. (45) 



Lastly, one can write (l44l) in normal form using the change of variables 6„ = 
/^n — ff/^n- The normal form of (jH]) at order 3 reads 

(3n+l - 2(3n + /5„-l = -^l Pn + B + h.O.t. (46) 

with 

5 = c, + | = ^(0^)(O)-^(V(^)(O))2). (47) 



This yields the explicit form fIzOl) of the reduced recurrence relation in the au- 
tonomous case {e} = {rj} = {7} = {k} = 0. 
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2.3.3 Homogeneous case near the upper phonon band edge 



In this section we consider the case ~ cui = v4+T]^, in the autonomous case when 
{e} = {rj} = {7} = {k} = 0. Equation ( ITTil takes the form ( 1271) . where L = A^^ is 
given by ([9]). The spectrum of L on the unit circle consists in a double non semi- 
simple eigenvalue —1, and the centre space is again spanned by = (cost,0), 
Vy = (O,cost). 

For fi = u'^ — uf in some neighbourhood A of 0, there exists a smooth two- 
dimensional local centre manifold A^^ C D locally invariant under L + N{.,fi), T, 
R and having the form (l29l) . The function ip having the centre manifold as its graph 
has the form fl30|) and shares the property fl3T|) . For /i ^ 0, the centre manifold 
Ai^ contains all solutions Yn of (I27p staying in a sufficiently small neighbourhood 
of y = in D for all n £ Z. Their coordinates (a„, 6„) on Ai^ are then given by 
a two-dimensional mapping, which determines all small amplitude solutions when 
0. 

The operator L has the following structure in the basis (V^, Vy) 

1 



and the reduced mapping is given by 



bn+1 J \ bn 



where 



6 ) ~ ( -a-26 + r(a,6,/i) ' ^^^^ 



and r is defined by ( IMI) . Since the reduced mapping inherits the symmetries of ( !27[) . 

commutes with Tjx^ = — / hence r(— a, — 6, /i) = — r(a, 6, /i). Moreover, fHHI) is 
reversible with respect to the symmetry R{a,b) = {b,a), which yields the identity 

r(a, b, /i) = r(— a — 26 -|- r(a, 6, /i), 6, /x). 

This implies r(a, 0, /i) = and 

r(a, 6, /i) = — 6/x + Ci6'^ + C2a6^ + -C2a% + h.o.t, (50) 

where the coefficients Ci, C2 have to be determined. 

For this purpose, we first compute the leading order terms in the Taylor expan- 
sion of ip at (a, 6, /i) = 0, using the fact that M.^ is locally invariant under L+iV(., /x). 
Equation (l36ll yields 

V9(-a - 26 + r(a, 6, /i), 6, yu) = (p{a, b, /x), (51) 
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(p{b, —a — 2b + r{a, b, fi), fj,) = 

{oof^ + 2 + fi2)<^(a, b, fi) - ^{b, a, fi) + {I - Pc) [yU^ + W] {b cost + (^(a, b, fi)). 

(52) 

The Taylor expansion of (p at order 2 takes the following form (due to the symmetry 
property (ISTl) ) 

ip{a, b, = (poiibfi + ^v^iioa^ + '■Piwab + ipo2ob^ + h.o.t. (53) 

By an identification procedure we now compute the coefficients (ppgr in ( l53l) . using 
( l52i) and the expansion ( HOl) . Identification at order bfi gives 

hence ipou = since y^oii ^ H^- Identification at order ab leads to 

1 

f020 = -{ujI— + 2 + fi2)(/?iio (54) 
and identification at order 6^ yields 

^110 + {^l^ - 2 + fi')<^020 = -^r^V^^'Ho) cos^t. (55) 



Reporting (IMI) in (1551) gives 



(^?^ + fiYv'no = -2fiV(3)(0) cos^t 



and consequently 

1 fi' 

(16 + 3^2) 



<^iio = -^^'H0)(— + _^^, cos(2t)), 



As a conclusion, we obtain 

y.(a,6,/i) = -r(3)(0)(^ + ^^^^cos(2t))(a6+ia2) 

-i^V(3)(0)( J, + A + (_|_, _ eos (2t)) 62 + h.o.t. 

(56) 

We now compute the two-dimensional mapping giving the coordinates (a„, 6„) 
of the solutions on Jvi^. Equation (jlHl) can be written 

ttn+i = bn, bn+1 + 26„ + bn^i = r(6„_i, 6„, yu). (57) 
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Using ([MD, dSHD and (HOD yields in equation (fSUD 

. Itf [V...(0) - (V->(0))^«'^(^ + ij + - ^)]. (58) 

e,._.-(^(3.(0)f[ijH-i^^^]. (59, 
The transformation 6„ = /5„ — y|/5^ yields the normal form of fl57l) of order 3 

Pn+l + Wn + = -/i + 5 + h.O.t. (60) 

with B defined by fl25p . This yields the explicit form (1231) of the reduced recurrence 
relation in the autonomous case {e} = {77} = {7} = {k} = 0. 

2.3.4 Inhomogeneous cases 

Using the normal form computations performed in the above sections for the au- 
tonomous case, one can obtain by perturbation the principal part (1201) (or (12^ ) of 
the normal form for uj ^ UJ2 (or uj ^ uJi) m. the nonautonomous case. In what follows 
this computation is described for uj ^ uj2i the treatment for u; ~ cji being completely 
similar. 

Theorem [3] is obtained by applying the reduction theorem [1] to the first order 
system f fTTj) . According to theorem [T]-i), small amplitude solutions F„ = {zn.Un) of 
( 1 11 1) have the following form for small {X} E E 

Yn = (fln, K) cost + ^(a„, bn, rn{\}), (61) 

where \&(a, b, {A}) = ip{{a, b) cost, {A}) G Dh and ip denotes the reduction function 
of theorem [H In the sequel we shall note = {"^i, "^2)- 

Let us compute the explicit form of the reduced map ffT^ . For this purpose, one 
has to use the ansatz flUTl) in equation (fTTl) and project the latter on the Fourier 
mode cost. Setting F„(a, 6, {A}) cost = /„((a,6) cost, {A}), the reduced map (fT4l) 
becomes 

bn+1 



(62) 



^"(■.W)(:) = (_„ + 2, (63) 

where (recall {A} = (/i, {e}, {77}, {7}, {k})) 

r,(a,6,{A}) = -{Q^n + + ^n)) b + + Vn){l + In) - l]b + Kn{b - a) 
+ P*W{b cost + ^2(a, b, r^{X}),Vn, 7n) 

(64) 
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and the function W is defined by f fT^ . 

Since /„(., {A}) commutes with T and T^Xc = the map {A}) commutes 
with —I. We have consequently 

r„(a,6,{A}) = -(l^^e^ + + e„)) 6 + + + _ 1]5 + ^.^(5 _ a) 
+ ci 6^ + C2 + C3 
+ 0(||(a,6)|n|{A}|U+||(a,6)f), 

(65) 

where the coefficients q need to be determined. Now, since rn[a,h, (yU, 0, 0, 0, 0)] = 
r{a,b,fj,) in the homogeneous case (see section [2.3.21) . we have C3 = — and ci, 
C2 are defined by ( 145|) . Consequently we have computed the principal part of the 
reduced equation (!62|) in the nonautonomous case. 

To obtain the normal form of of order three we now define P{f3) = (3 — ^(3^ 
and consider as in section 12.3.21 

an = P{an), hn = P{Pn)- 

This yields the normal form of fl62p of order 3 given in equation fl2UI) . 
Moreover, the small amplitude solutions of ([6]) have the form 

yn={Pn-^Pl) COSt + Vl>2(P(/5„_i),P(/?„),r„{A}), 

therefore the reduction function of theorem[3]is given by 0(a, /?, {A}) = — ff/?^ cost+ 
\E'2(-P(«), P{P), {^})- Note that the reduction function has a component along the 
Fourier mode cost after the normal form transformation. 

3 Exact periodic solutions for an homogeneous 
lattice 

Here we consider the case of the homogeneous Klein-Gordon lattice ([3]), which leads 
us to system ([6]) with {e} = {rj} = {7} = {k} = 0. Breather solutions have been 
proved to exist by MacKay and Aubry |MA94] for system ([3]) with small values of 
the coupling parameter k and nonresonant breather frequencies. Here we prove the 
existence of small amplitude breathers for arbitrary values of k in some cases and 
frequencies close to the phonon band edges (see theorems O-i) , E]-i) and [7] below). 
We also prove the existence of dark breather solutions, which converge towards a 
nonlinear standing wave as n ^ ±00 and have a much smaller amplitude at the 
centre of the chain. 

Let us start with the case cu ~ O and fl > 2/ \/3 in (Q. By theorem [3l small am- 
plitude solutions of (EI) in Hi are determined by the recurrence relation ( |20l) . This 
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recurrence becomes autonomous for an homogeneous lattice and takes the form 
It is important to note that the invariance n —>■ —n of ^ in the homogeneous case is 
inherited by ( H6|) (see |Jam03] . section 5.2 and theorem 2). This invariance imphes 
that the two-dimensional map i— >• (/?„,/3„+i) is reversible. Bifurcations of 
small amplitude homoclinic and heteroclinic solutions have been studied in |Jam03j 
(section 6.2.3) for this class of maps. This yields the following result for the recur- 
rence relation 



Lemma 2 Assume n > 2/V3 andB = f{ V^^^^O) - |(V(3)(o))2 ) q. For fi ^ 0, 
the recurrence relation [J^ j has the following solutions. 

i) For fi < and B < 0, (^) has at least two homoclinic solutions (3\, (3^ (and 



also —Pn' ~Pn) that \imn^±ooPn = 0- These solutions have the symmetries 
P\+, = PI P\ = PI and satisfy < /5; < C a^^'^^ with a, = 1 + 0(|/i|i/2) > 
1. 

ii) If fi and B have the same sign, has two symmetric fixed points ±/5* = 

0(|/.|V2). 

Hi) For fi > and B > 0, [J^ j has at least two heteroclinic solutions P^ (and 
also —Pn, ~Pn) ^^^^ ^^^^ lim„^±oo/?n = Thcsc solutions have the symmetries 

P^nJ^i = —Pn and = —Pn- Morcovcr, are 0(/i^/^) as n ^ ±oo and 0{ix) 

for bounded values ofn. 

Note that for B > and fi < {fi ^ 0), pUj) has no small amplitude homoclinic 
solution to 0. 

For fi < 0, typical plots of the stable and unstable manifolds of the fixed point 
Pn = are shown in figure [31 page [HUl (nonintersecting case B > 0) and in figures [7] 
and fT2] pages 1391 and 1531 (intersecting case B <0). 

Theorem [H] ensures that each solution in lemma [2] corresponds to a solution 
Vn of ® given by 

yi,{t) = Pi^ cost + 0(/5;_i, PI {fi, 0, 0, 0, 0)) (66) 

with tu^ = fi^ + yU in (ini). The following result follows (the symmetries of y^ are due 
to the symmetries of P^ described in lemma E]). 

Theorem 5 Fix {e} = {77} = {7} = {n} = in equation Assume Q > 2/V3 
and b = V^'^\0) - |(l^(^)(0))2 ^ 0. For 00 ^ n, problem ^ has the following 
solutions with yn G for all n eTj. 

i) For oj < VL and b < ^ has at least two homoclinic solutions y\, y^ (and 
also y\{t + tt), y^(t + 71)) such that lim„^-|-oo llz/nllHl = 0- These solutions satisfy 
y-n+i = yi, y-n = yl the form 

= cost + 0(1^ -r^l), (67) 



28 



where < Pi < C \uj - n\^/^ and ai = 1 + 0{\uj - > 1. Solutions yl, yl 

correspond to small amplitude breathers with a slow exponential decay as n ^ ±oo. 
ii) Ifuj — Q and b have the same sign, ^ admits a solution y^ G independent of 
n, corresponding to collective in-phase oscillations. It has the formy^{t) = (3* cost+ 
0{\u; - n\) and f3* = 0{\uj - 

Hi) For LJ > Q and b > 0, ^ has at least two heteroclinic solutions y^, y^ 

(and also yl^{t + vr), y^{t + tt) j such that lim„^_oo — + tt) ||^2 = and 

# 

lim„_^+oo WVn ~ y^Wu^ = 0- These solutions satisfy y^^_^_i(t) = ^^(t+vr) andytnit) = 

y^{t+7r). Moreover, their norms ||?/^||jy2, ||?/^||j|i-2 areO{{uj — ^iy^'^) as n ^ ±oo and 

0{{uj — VL)) for bounded values of n. Solutions y^, correspond to small amplitude 
dark breathers. 

In addition, note that for 6 > there exists no small amplitude discrete breather 
yn G with u < fl and u ^ fl (since ( HUj) has no small amplitude solution 
homoclinic to 0). 

Now we consider the case to ^ with tOc = v^4 + ^7^. In that case, equation 
(in]) can be locally reduced to the recurrence relation fpH|) . which becomes again 
autonomous if {e} = {77} = {7} = {k} = and has the invariance n —n. This 
class of recurrence relations has been studied in |Jam03] (section 6.2.3, lemma 7) 
to which we refer for details. In addition one can note that the recurrence ( I6OI) can 
be recast in the form (1461) by setting = (—1 The following result for the 

recurrence relation flUUj) follows. 

Lemma 3 Assume B = V^^\0) + {V^^\0)f{j^^ - 2) ) ^ 0. For 12 ^ 0, the 
recurrence relation IfUI^) has the following solutions. 

i) For fi > and B > 0, has at least two homoclinic solutions (3\, (3"^ (and 
also —Pn) that lim„^-|-oo = 0. These solutions have the symmetries 
P-n+i = -PL P-n = Pi and satisfy < (-!)"/?; < C/i^/^ ^^^^ |^^| ^ 
l + 0(|/i|V2) > 1. 

ii) If fi and B have the same sign, l[U(]\) has a period 2 solution /3° = (— with 

Hi) For fi < and B < 0, f6U\} has at least two heteroclinic solutions P^, P^ (and also 
—Pn, —Pn) such that lim„^±oo \Pn T Pn\ = 0- These solutions have the symmetries 
P^n+i = Pi P-n = ~Pn- Morcovcr, P^, P^ are 0(|/i|^/^) as n ^ ±00 and 0(|/i|) 
for bounded values ofn. 

In addition, for < and fi > {fi ^ 0) problem (I6UI) has no small amplitude 
homoclinic solution to 0. As above, the solutions of the reduced recurrence relation 
provided by lemma [3] yield the following solutions of ([6]). 
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Theorem 6 Fix {e} = {r/} = {7} = {k} = in equation Assume b = 

V^^\0) + {V^'^\0)f{j^^ - 2) ^ 0. For uj ^ 00^ = V^+W, problem ^ has the 
following solutions with Un € H"^ for all n eTj. 

i) For ijj > ijJc and b > 0, ^ has at least two homoclinic solutions y\, (and 
also y\(t + Tx), y'^it + rr)) such that lim„^-|-oo ||2/nll//2 = 0. These solutions satisfy 
V-n+iit) =yni't + TT), y\ = yl and have the form 

y^ = 131, cost + 0{\uj-uj,\), (68) 

where Q < (-1)"/?; < C (cj - cuj^/^ |^^|-|n| |^^| ^ i ^ o{{ijj - uj^f'^) > 1. 
Solutions y\, y^ correspond to small amplitude breathers with a slow exponential 
decay as n ^ ±00. 

ii) If uj — uJc and b have the same sign, ^ admits a solution y^ being 2-periodic in 
n, corresponding to out-of-phase oscillations. It has the form y^{t) = y{t + mT) with 
y{t) = (3* cost + 0{\uj - uj,\) (y G Hi) and /3* = 0{\uj - ujc\^/^). 

Hi) For UJ < ijjf. and 6 < 0, ([^ has at least two heteroclinic solutions y\, y^ 
(and also y^(t + n), y!^(t + n)) such that lim„^_oo — + vr)||^2 = and 

lim„_,+oo llz/n — y^ll/fi = 0. These solutions satisfy y'Ln+i = Un '^"■'^ V-ni't) = 
y^{t+n). Moreover, their norms \\yn\\H^ , Ill/nil//^ are 0{\uj — uJc\^/'^) as n ^ ±00 and 
0{\uj — ujc\) for bounded values ofn. Solutions y^, y^ correspond to small amplitude 
dark breathers. 

In addition, for 6 < there exists no small amplitude discrete breather y^ € 
with UJ > UJc and u ^ uJc- 

It is worthwhile mentioning that approximate breather solutions of (j6]) can be 
obtained in the form of modulated plane waves, using multiscale expansions (see 
|GM04j and references therein), where the error can be controlled over finite time 
intervals. The envelope of a modulated wave satisfies the nonlinear Schrodinger 
(NLS) equation, and does not propagate along the chain when a plane wave with 
wavenumber g = or g = tt is modulated (its group velocity vanishes). In these 
two cases the NLS equation is focusing (i.e. time-periodic and spatially localized 
solutions exist) when 6 < and 6 > respectively, which coincides (according to 
theorems and [H]) with the parameter values for which exact breathers exist. 

In addition, as shown in reference |Fla96] the condition 6 < leads to the 
instability of nonlinear standing waves with wavenumber g = 0. If periodic boundary 
conditions are considered, these standing waves become unstable above a critical 
energy via a tangent bifurcation. When the lattice period tends to infinity, the 
energy threshold goes to and bifurcating solutions are slowly spatially modulated. 
The same result has been obtained for standing waves with q = tt when V is even 
and 6 > 0. 
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In what follows we reformulate the results with respect to the unsealed original 
system For conciseness we only describe breather bifurcations, but conditions 
for dark breather bifurcations are easily deduced from theorems [5] and El We express 
the condition 6 > of theorem [6] in a different way using the relation 

16 (V(3)(0))^ 
3 16 + 3^]2• 

In addition, as the rescaled potential y of (E]) is replaced by the original potential 
^ of ([3]), coefficients b and b are simply replaced hj h = a~^b and h = a~%. 

Theorem 7 Consider the Klein-Gordon lattice where the on-site potential V 
satzsfiesV'iO) = 0, V"{0) = 1 andm,d,a,k> 0. Assume h = \/W(0)-|(\/(3)(0))2 ^ 
and note fi^ = a^d/k, cu^.„ = a'^d/m, cj^„^. = (a'^d + 4:k)/m and H = 1/^(0) - 

i) If h < and Q"^ > 4/3, system ^ admits two families of breather solutions x\, x,^ 
parametrized by their frequency uo (in addition to phase shift), where uo ^ ujmin ond 
uj < uJmin- These solutions satisfy x\^j^i = x\ and x^^ = x"^ and decay exponentially 
as n ^ ±oo. As UJ uJmin, the amplitude of oscillations and the exponential rate 
of decay are 0{\uj — uJminl^^'^)- The breather profile is a slow modulation of a linear 
mode with wavenumber q = 0. 

a) If h > and Q"^ > —16H/{3h), system ^ admits two families of breather 
solutions x\, x\ parametrized by their frequency u) (in addition to phase shift), where 
ijj ^ ijJmax (ind UJ > ujmax- Thcsc solutions satisfy x^„_^_-^(t) = x\(t -\- tt / lj) , x'^_^ = x^ 
and decay exponentially as n ^ ±oo. As u ^ uJmax, the amplitude of oscillations 
and the exponential rate of decay are 0{\u — uJmax\^^'^)- The breather profile is a slow 
modulation of a linear mode with wavenumber q = tt. 

To interpret the conditions on the on-site potential V in properties i) and ii), it 
is interesting to note that V is soft for h < and hard for h > near the origin (i.e. 
the period of small oscillations in this potential respectively increases or decreases 
with amplitude). The condition on Q in property i) corresponds to a nonresonance 
condition, i.e. it ensures that no multiple of lj lies in the phonon band [uJmin, ^max] 
for UJ ^ UJmin- The condition on Q in property ii) is of different nature and is 
equivalent (with the condition /i > 0) to fixing h > 0. 

Discrete breathers were known to exist in Klein-Gordon lattices for small cou- 
pling k after the work of MacKay and Aubry |MA94j . Theorem [7] considerably 
enlarges the domain of breather existence, with the limitation that it only provides 
small amplitude solutions. In particular, it is interesting to note that small ampli- 
tude breathers of property ii) exist for all values of k ii H > 0. 
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4 Normal form analysis for inhomogeneous lat- 
tices 



According to theorem [3l small amplitude solutions of (Q in are described (for 
small inhomogeneities and frequencies close to the phonon band edges) by finite- 
dimensional nonautonomous recurrence relations. In what follows we only consider 
the case u ^ fl, the situation when u ^ a/4 + fi^ yielding similar phenomena. At 
leading order, the reduced recurrence relation (l20l) can be approximated by 



(69) 

Different kinds of techniques can be employed to obtain homoclinic solutions of flU^ . 
One can use variational methods for asymptotically periodic sequences |PZ01j (see 
also [Wei99] in the homogeneous case), or proceed by perturbation near an uncou- 
pled limit (also denoted anti-continuous or anti-integrable limit) where Q and B are 
large (see |AA90l IAub95l IABK04j and section 9 of jMA94] ). Existence results of 



localized solutions are also available in |AF88j for disordered defect sequences. An- 
other approach is to start from a known uniformly hyperbolic homoclinic solution 
in the homogeneous case, which persists for small inhomogeneities by the implicit 
function theorem, and obtain estimates for defect sizes allowing persistence (see 
the technique developed by Bishnani and MacKay |BM03j ). Interesting related re- 
sults on the structural stability of discrete dynamical systems under nonautonomous 
perturbations can be found in reference |Fra74j . 

With a different point of view, we develop here a dynamical system technique, 
valid for a finite number of defects, which allows to analyze bifurcations of homoclinic 
solutions as defects are varied (see sections Wl] and Hl2l) . For an isolated defect we 
put in evidence, near critical defect values, bifurcations of new homoclinic solutions 
(having no counterpart in the homogeneous system) or disappearance of homoclinic 
solutions existing in the homogeneous case. Our method is also generalized to a 
finite number of defects, with the counterpart that ( |69l) is modified by suitable 
higher order terms depending on the defect sequence (this procedure only provides 
approximate solutions of (|69|) ). However this does not constitute a strong limitation 
since the full reduced equation fl^Uj) is itself a higher order perturbation of flU^ . Note 
that equation flUUj) is valid (according to theorem [3]) for small defect sizes and /i ~ 0, 
where the parameter /i determines for /i < the (weak) degree of hyperbolicity of 
the fixed point = in the homogeneous case. Our analysis does not impose 
conditions on the relative sizes of these parameters. 

In order to obtain exact breather solutions of (jH]) via theorem [31 it would be 
necessary to proceed in two steps. The first step is the one described above, where 
exact or approximate homoclinic solutions are obtained for the truncated problem 
( l69l) . The second step is to show that these solutions persist for the complete 
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equation fl^ as higher order terms are added. For this purpose a typical procedure 
would be to solve (1201) using the contraction mapping theorem in the neighbourhood 
of an exact or approximate solution of ( !69l) . 

In this paper we shall not examine the persistence of solutions for the complete 
equation ( l20l) . Instead we shall later compare approximate solutions yn(t) ~ Pn cost 
(deduced from fl69p and theorem [S]) to numerically computed solutions of the original 
problem (see section E]). This will allow us to study the validity of approximation 
(1691) far from the small amplitude limit and as inhomogeneities become larger. 

Note that other interesting bifurcations can exist when impurities act at a purely 
nonlinear level (see [SKRCOll IKKK03 ] for some examples in spatially discrete or con- 
tinuous systems). This would correspond to the situation when the on-site potential 
in has an harmonic part independent on n, whereas higher order terms are in- 
homogeneous. The subsequent analysis of the reduced recurrence relation would 
be quite different, and in particular the method developed in sections l4.ll and IT2] 
(based on a linear deformation of the unstable manifold) would not apply. 

4.1 Case of a single mass defect 

We start with the simplest case when the coefficients of ( l69l) are constant, except at 
n = where their value changes. To fix the idea we consider the case of a single mass 
defect in equation (0), i.e. Dn = d, = k, = a, Mn = 1 + m„ = mo^no- 
The case when all lattice parameters are allowed to vary over a finite number of 
sites will be considered in section 14.21 

For equation ^ the above assumption yields ?7n = 7n = = and e„ = mo^no- 
Equation (lUOj) reads (recall oj"^ = Vl? + ji) 



Pn+l - Wn + /?n-l = -(o^^mQ^nO + /i) /?n + ^ 

Setting = an and f/„ = (a„,/3„)"^, equation (!70l) can be rewritten 



(70) 




(71) 



where 




an + 2/3„ + (02 _ uj^)i3^ + BPl 




(72) 



One has in particular 



Ui = A{Lj,mo)G^{Uo) 



(73) 



where the linear transformation 




(74) 
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corresponds to a linear shear. Note that the axis a = consists of fixed points of 
A{uj,mo). 

It is worthwhile to notice that the map is reversible under the symmetry 
R : (a,/?) I— i> (/?,«), i.e. G^^ o R = RG~^. In other words, if Un is a solution of 
( 17Ti) for mo = then R U-n is also solution. This property is due to the fact that 
equation (ITOjl with mo = has the invariances n — >■ n + 1 and n —n. Obviously 
the latter invariance still exists for mg 7^ 0. Consequently, for all mo € M, if t/„ is a 
solution of ( ffB then R U-n+i is also solution. 

Now we shall use a geometrical argument to find homoclinic orbits to for 
equation ( ITOll . In the sequel we consider the stable manifold iy*(0) of the fixed 
point {a, (3) = of G^, and its unstable manifold iy"(0), both existing for u < Q. 
The following result follows immediately. 



Lemma 4 For < u < Q, equation [7W possesses an homoclinic orbit to if and 



onlyifW'^{0) and A{uj,mQ)W^{0) intersect. 

In addition it is useful to notice that Vr'*(0) and H^"(0) are exchanged by the 
reversibility symmetry R. 



4.1.1 Linear case 

As a simple illustration, consider the linear case when V is harmonic, in which 
B = 0. Equation (ITUj) reads 

Pn+i - 2/5„ + Pn-1 = -(u;^mo(5„o + /i) (75) 

In that case, Vr*(0) and W'^{0) correspond respectively to the stable and the unstable 
eigenspace of a linear mapping in M^. The situation is sketched in figure H] below. 
The corresponding stable eigenvalue a G (0, 1) is given by 

^ = l-^-^(/i^-4/i)^^ f^ = u;'-n'<0, (76) 

and W^{0) is the line /3 = era, W'^{0) being the line (3 = a~^a. For fixed u < ^l, 
A{uj, mo) maps the unstable eigenspace on the stable eigenspace if and only if mo > 
(mass is increased at the defect) and 

mo = mi{uj), (77) 

where 

mi{u;) = ^{a-' -a). (78) 
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Figure 4: Linear case {B = 0). (Left panel) Stable manifold (in the half plane 
a > 0) and images of the unstable manifold by A{uj,mQ) for mo = 0.005, tuq = mi 
and mo = m/ + 0.005. (Right panel) Homoclinic orbit to for mo = mi. In both 
panels we have fixed Q = 10, u = 9.99, which implies m/ = 0.0092. 



Now keeping fixed mo > 0, condition 0771) can be rewritten u = C(j;(mo), where 
(for uj <^l) 

cof = -^[n^ + 2 - (4 + m2fi2(f]2 ^ 4))i/2 ]^ moj^l, 
1 mQ 

The solutions of fl75|) homoclinic to are spanned by /5„ = a''^', and the correspond- 
ing solutions of in the linear case read yn{t) = /3„cost with u = ui{mQ). One 
recovers a classical result, i.e. if mass is increased at the defect then the linear 
localized mode frequency lies below the phonon band and its frequency is given by 

UJl. 

Now let us consider the effects of nonlinear terms. For this purpose we start with 
the simplest case of a hard potential, i.e. B > 0. The situation when 5 < is far 
more complex and will be investigated later. 



4.1.2 Nonlinear defect modes for hard on-site potentials 

If i? > 0, iy*(0) and H^" (0) do not intersect (except at the origin) ior < cu < Vl. 
Indeed, one can show by induction that \Pn\ > \Pn-i\ > for any nontrivial orbit 
on Vr"(0), which implies VT" (0) lies inside the sector formed by the lines a = (3 and 
a = 0. In the same way, W^''(0) lies inside the sector formed by the lines a = (3 
and (3 = hence it does not intersect iy"(0). The above property also implies that 
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Vr"(0) can be defined (globally) as the graph a = g{(3) of an increasing function g, 
and the same holds true for Ty*(0) = RW^{Q) on which (3 = g{a). 

For fixed u G (0,f2), the local unstable manifold can be approximated by 
a = g{p) = aP + bp^ + 0{\p\^), with b = a\a^ - a~^)~^B < (this coefficient 
can be computed by a classical identification procedure, using the fact that W^{0) 
is invariant under G^^). Consequently, 1^(0) and Vr"(0) have the local shape rep- 
resented in figure The same situation occurs in the limit u ^ Q (one can locally 
approximate the map G^^ up to any order in U, fi using the time-one map of an 
integrable fiow |AP90j . which allows to determine the shape of iy*(0) and iy"(0) 
close to Un = 0). 

In the case when mo < 0, the curves 1^(0) and A{uj,'mQ)W'^{0) do not intersect 
(A(a;, mo)H^"(0) remains inside the sector formed by the lines a = P and a = 0). 
However, W^{0) and A{uj,mQ)W^{0) intersect if uiq > provided 

mo > mi{u), (79) 

which means there exists an orbit homoclinic to for equation fl70l) . This property 
is clear for mo ~ rrii where there exists a unique intersection point (in the half plane 
a > 0) close to f/ = due to the local shape of iy*(0) and W^{0). Moreover, we 
numerically find a unique intersection point for all values of mo satisfying (!79|) . 
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Figure 5: Case -B > and uj < VL. Stable and unstable manifolds (in the half 
plane a > 0), and image of the unstable manifold by ^4(^7, mo) for mo = m^ and 
mo = 0.05 > mi. We have fixed = 10 and uj = 9.9. 

Condition (1791) is equivalent to 

uJi<u <Q. (80) 

The amplitude of the homoclinic orbit is 0{^yuJ — uji) as mo is fixed and uj ^ ui 
(at the limit y4(u;, mo)W^"(0) and 1^(0) become tangent at the origin), and its 



36 



spatial decay rate a is given by equation flTUl) . This homoclinic orbit corresponds 
for equation ([H]) to a nonlinear analogue of the above mentioned linear localized 
mode. This solution can be approximated by Unif) ~ PnCost for mo ~ and 
the frequency uj varies with amplitude contrarily to the linear case. Note that the 
existence of this nonlinear defect mode can be obtained in a standard way from 
equation in the small amplitude limit, using the Lyapunov centre theorem (in its 
infinite-dimensional version). The Lyapunov family of periodic orbits exists when no 
multiple of ui lies in the phonon band [fi, (fi^ + 4)^^^], which is true in our parameter 
regime where Q, > 2/ y/S and mo ~ . 

Lastly, let us notice that the above homoclinic orbit possesses the symmetry 
P-n = Pn, or equivalently i?f/_„+i = Un- It suffices to check the latter relation for 
n = to prove it for any n, since both solutions RU^n+i and f/„ coincide if they 
satisfy the same initial contidion. Since Uq lies on the unstable manifold we have 
ao = g{Po), and in the same way /?i = g{ai) since Ui lies on the stable manifold. 
Since by definition ai = Po, this implies ao = Pi and thus RUi = Uq- Using the 
properties Ui = G^{Uo) — u'^rriQ^O, Pq)^ and RUi = Uq we also deduce the relations 

2ao = [2 + ~ uj\mo + l)]po + B P^ (81) 
2pi = [2 + - uj\mo + l)]ai + B al (82) 

which are useful in particular for the numerical computation of Uq, Ui. 
4.1.3 Nonlinear defect mode with algebraic decay 

In the situation of section 14.1.21 {B > 0), the case when mo is fixed and uj Q 
deserves a special attention. Indeed, the homoclinic orbit (a„,/5„) converges in this 
limit towards a solution having an algebraic decay as n ^ ±oo. 

More precisely, if u = Q and mo > then iy*(0) and A{uj,mo)W^{0) intersect 
at a unique point (ai. Pi) in the half plane a > (see figure E]). This can be checked 
analytically for mo ~ and {a,P) ~ 0, since the unstable manifold can be locally 
parametrized by a = /? — (-8/2)^/^/5^ + 0(1/5^) in the half plane a > (this expansion 
follows from a classical identification procedure). Using this relation for {ao, Pq) in 
conjunction with (1811) we find as mo — ^ 

Po = n\2B)-^/^mo + Oiml). (83) 

Note that in this non-hyperbolic case, the function g having the unstable manifold 
as its graph is not at /5 = (in the half plane a < 0, one has a = /3 + (5/2)^/^/5^ + 
0(|/?|^) on the local unstable manifold). Far from the small amplitude limit, we have 
also checked numerically the existence and uniqueness of iy*(0) fl A{uj,mo)W^{0) 
in the half plane a > 0. 
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Figure 6: Case -B > and uj = VL. Stable and unstable manifolds (in the half plane 
a > 0), and image of the unstable manifold by A{uj,mQ) for mo = 0.05. We have 
fixed n = 10 in this example. 



Consequently, there exists a solution of ( !70l) homoclinic to for = f2 and 
mo > 0. This solution has an algebraic decay due to the fact that the origin is not 
any more an hyperbolic fixed point for uj = Vt. One can approximate the solution 
profile for mo ~ 0, using the fact that (170|) admits at both sides of n = a continuum 
limit. Indeed, setting 

(3n ~ mo/?(x), X = rriQn, (84) 
one obtains the following differential equation 

-4 = S/?^ X e (-00, 0) or (0, +oc), 
from which we deduce (multiply by (3' and integrate) 

^ = -sign(x) (5/2)1/2^2 (g5) 

since (3{x) ^ as x ^ ±oo. Using and (IH^ one obtains the following approx- 
imation of the homoclinic solution for mg ~ 

lY 2 , 

Pn ^ mo d — {mo\n\ + —) . (86) 

This yields an approximate solution yn{t) ~ /5„cost of (EI), corresponding to 
a breather with an algebraic decay and a frequency uj = Q a.t the bottom of the 
phonon band. 
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4.1.4 Case of soft on-site potentials 

Now let us make some considerations on the case B < (soft on-site potential V) 
which is far more complex. 

For mo = 0, < and B < 0, equation flTOj) possesses homoclinic solutions to 
0. This case has been analyzed in several references with different viewpoints and 
for different parameter ranges, see e.g. |HT99[ IABKn4l [PZOTj IWei99t I.TamOaj . 

The dynamics of the map Go; is rather complex due to the fact that the stable 
and unstable manifolds of the origin intersect transversally in general (see figures [7] 
and [12]). This implies the existence of an invariant Cantor set on which some iterate 

is topologically conjugate to a full shift on symbols |AP90] . which yields a 
rich variety of solutions and in particular an infinity of homoclinic orbits to 0. 

Among these different homoclinic orbits one can point out two particular ones 

= (a^,/?*)"^ {i = 1,2), corresponding for the Klein-Gordon chain to breathers 
solutions with a single hump near n = (site-centered or bond-centered). These 
solutions have been described in lemma[2]and theoremOin the small amplitude limit. 
The corresponding homoclinics U^, are reversible, i.e. they satisfy Rlf^n+i = 
[j3\ = (31) and i?f/i„+2 = iP-n+i = Pi)- In figure [3 the point with label 2 lying 
on the axis a = (3 corresponds to Ul, and the points with labels 3, 1 correspond to 
Uq,U1 respectively. Obviously any translation of f/^ generates a breather solution 
of ^ having its maximal amplitude near a different site. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 



Figure 7: First intersection points between the stable and unstable manifolds for 
parameters u = 9.9 (/x = —1.99) and B = —75. The dashed line depicts the image 
of the unstable manifold by the linear shear A{uj,mQ) for rriQ = 0.005. 

Now let us consider the situation when u is kept fixed and a small mass defect 
mo is introduced in (!70i) . As illustrated in figure [71 each of the above solutions is 
structurally stable. 
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For example, let us consider in figure [7] the intersection points 1,2,3 between 
iy"(0) and iy*(0). Each of these intersections persists (points l',2', 3' in figure [7]) 
as the linear shear A{uj,mQ) is applied to W^"(0) for mo ~ (dashed line in figure 
[7j). Let us examine the corresponding solutions of ( 17T1) and the related breather 
solutions of the Klein-Gordon model. 

We denote by = (a^, P^)'^ the solution of (!7T|) with initial data Uf at the point 
1'. This solution is homoclinic to according to lemma HI Repeating an argument 
of section H:.1.2l one can show that RU'^n+i = i-^- P-n = Pn- Consequently, 

corresponds to an (approximate) breather solution of ([6]) centered at the defect 
site n = 0. This solution is a small deformation of the site-centered breather y'^ of 
theorem O 

Now let us denote by the homoclinic solution of (1711) with initial data Ul 
at the point 2'. It corresponds to an (approximate) breather solution of ([6]), whose 
profile is a small deformation of the breather centered between n = and n = 1 
(see theorem E]). Since Ul does not belong to the line a = P {it lies at a distance 
0(|mo|)), the corresponding breather solution is not symmetric any more, which 
was expected since the atomic masses at n = 0, 1 are different. 

Lastly we note = (aj^,/3^)"^ the homoclinic solution of ( 17T1) with initial data 
Uf at the point 3'. Since Uf is O(|mo|)-close to Uq (point with label 3), is a 
small deformation of the solution U!^_i existing for tuq = 0. In other words. Iff 
corresponds to a small deformation of the breather y'^_i centered at n = 1. The 
mass defect at = breaks the mirror symmetry of the solution, since its amplitude 
has only the imperfect symmetry P'^n+i ~ Pn+i = C^d^^ol) for n ^ 0. 

A more delicate question concerns the continuation and the possible bifurcations 
of the above homoclinic solutions as mo is further varied. The evolution of U^, f/^, 
depends on the structure of the homoclinic windings near Ul,Uf,U^. Numerically 
we find that the lobes formed near these points by the stable and unstable manifolds 
have the structure shown in figure [3 These manifolds windings can be analytically 
approximated as explained in reference [HT99j (section 3.5) or [HRGB96] (section 
4). 

At a critical value mo = mc{uj) > 0, the points with label 2' and 3' on iy*(0) fl 
A{uj, mo)Vr"(0) colhde as W^*(0) and A{u!, mo)W^{0) become tangent. Consequently 
the solutions and disappear through a tangent bifurcation above this critical 
value of mo. Obviously, since we consider the truncated map flTTj) instead of the full 
recurrence relation (|20l) . these solutions only correspond to approximate breather 
solutions of (E]). However we shall check numerically that the corresponding tangent 
bifurcation occurs for nearby breather solutions of ([6]), at a critical value close to 
mc{w) (see section |5]). 

In what follows we give a simple method to estimate mduj), which is based 
on a simple approximation of W^{0). Let us consider a cubic approximation 
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of the local unstable manifold of figure U\ parametrized hj /3 = Xa — (? . The 
coefficient c depends on /i and B and need not be specified in what follows (a value 
of c suitable when A is large is computed in |HRGB96] . equation (60)). We note 
A = o^^ = 1 — /i/2 + A//i^ — 4yu/2 the unstable eigenvalue. We have 

/3 = Aoa - 0? (87) 

on the curve A{u)^ '^^^appi where A— cu^mo = Aq. By symmetry we can approximate 
the local stable manifold using the curve parametrized by 

a = A/3-c2/?3. (88) 
The curves A{u)^raQ)W'^^^ and W^^^ become tangent at (a,/?) when in addition 

{\-'ic^f){\o-'ic^o?) = \. (89) 

In order to compute mo = vnc as a function of tu, or, equivalently, the correspond- 
ing value of Ao as a function of A, one has to solve the nonlinear system (!871) -( l88l) - (!89l) 
with respect to a, /5, Aq, which yields a solution depending on A. Instead of using A 
it is practical to parametrize the solutions by t = This yields 

a = ^{t^-S^I^ /3=^ ft +1)1/2 



3 1 3 1 3 

^°-2^+2t^' ^-2t + 2'^- 
Since /x = 2 — A — A~i and mo = (A — Aq) + /i)~i it follows 



3 + t^ 2t 

u = 2 -, 90 

mo = ^(t-^)^(fi^ + /i)-^ (91) 

Given a value of yU G (— —1/2), one can approximate rric by the value of mo given 
by equations CT-CT. 

For example, in the case numerically studied in figure [7| we have u; = 9.9 and 
^ = -1.99. Consequently A ^ 3.721, t ^ 1.7935 and Aq ^ 2.777, which yields 
rric ~ 0.009632. A numerical study of the map yields rric G (0.00963,0.00964), and 
consequently our approximation works very well in this parameter regime. Moreover, 
the approximation is extremely close to the actual value of mo at which a tangent 
bifurcation occurs between the corresponding breather solutions of the Klein-Gordon 
system (numerically we again find mo € (0.00963,0.00964), see section [5] for more 
details). 
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Despite it gives precise numerical results in a certain parameter range, the ap- 
proximation flUUl) - fl^JTl) is not always valid. Indeed, the parameter regime fi > —1/2 
is not described within this approximation. Moreover, one can check that VF^^ in- 
tersects W^pp on the line a = P with an angle depending solely on A, and not on 
the coefficient B (in particular, W^pp and W^pp become tangent for A = 2). This 
problem could be solved by adding a quintic term d ot" in equation (IH7I) . 

The intersection point with label 1' between iy*(0) and y4(ti;, mo)Vr"(0) persists 
for < mo < raiiuS)^ or equivalently < < ci;;(mo), and consequently the reversible 
homoclinic solution [/^ exists within this parameter range. At = cj^, this solution 
disappears through a pitchfork bifurcation with the symmetric solution — t/^ (the 
amplitude of the homoclinic orbit is 0{^\/io — ui) as mo is fixed and uj —>■ uji). This 
homoclinic orbit corresponds for equation ([6]) to a nonlinear analogue of the linear 
localized mode of section I4.1.1[ As noticed in section 14.1.21 the existence of small 
amplitude exact solutions of this type (with oj ^ oji) can be obtained in a classical 
way using the Lyapunov centre theorem, when uji does not resonate with the phonon 
spectrum. For mo ~ and u; ~ cu/, the breather solution of ^ can be approximated 
by Unif) ~ Pn cost, and the frequency u varies with amplitude and lies below ui. 

More generally, the evolution of the set y4(a;, mo)W^"(0) fl iy*(0) as mo varies is 
very complex, due to the complex shape of the stable and unstable manifolds and 
the complicated structure of their intersection (see figure [T2i) . In section [5] we shall 
give some additional examples of breather bifurcations which can be deduced from 
the fine structure of the stable and unstable manifolds. 

Note that previous studies have examined, for certain families of reversible two- 
dimensional maps, how parameter changes modify the intersections between the 
stable and unstable manifolds of the origin and the associated set of homoclinic so- 
lutions |BCKRBW00l m^TCMOGj . These (autonomous) maps are directly obtained 
from the discrete nonlinear Schrodinger equation or generalized versions (due to 
their phase invariance), as one looks for oscillatory solutions with a single Fourier 
component. Althought we obtain similar types of tangent bifurcations as defects 
strength are varied, our situation is quite different since we are concerned with a 
nonautonomous map, where the impurity leads to consider a linear shear of the 
unstable manifold. 

4.2 Case of finitely many defects 

This section generalizes the analysis of the above section to the case when equation 
(E]) admits a finite number of inhomogeneities. More precisely, we assume in equation 
(E]) = 7n = = Cn = if \n\ > Tio + 1, for a given integer no > 0. Note that this 
assumption allows one to cover the case of an odd number of defects as well as an 
even number. 



42 



The situation is more complex than in section 14.1^ because studying homochnic 
solutions of flU^ leads to finding the intersections of the stable manifold with the 
image of the unstable manifold under a nonlinear transformation (see lemma [5] 
below). However, one can recover the linear case if one replaces the relevant spatial 
map by a suitable one, both being equal at leading order only. In what follows we 
shall develop this leading order theory, considering as higher order terms all terms 
being o{\\{an, Pn)f)- 

Equation ( l69l) reads 

/?„+! - 2(3n + Pn-l = {On - f^)f3n - ^^nPn-l + B (92) 

where On = ^"^{rjn + rjnln + In) — ^"^^n + «^n, oj"^ = Vt^ + fi. In the sequel we shall 
note e = \\ {0} + || {k} ||^^(z). 

Setting Pn-i = «n and Un = {an,Pn)'^, equation (lU^ can be rewritten 

Un+l = F^iUn), (93) 
^ ( -(1 + Kn)a +{/+0n~ + 5/33 ) • ^^^^ 

Noting F = for simplicity (see definition (17^ ). one can observe that 

F„ = (J + T„)F + 0(|«:„||/5|3), (95) 



T 





On + {fi- 2) 

^n ^n 

Note that higher order terms are absent from equation (I95p if k„ = 0. 
Since F„ = F for \n\ > uq + 1 one has the following property. 

Lemma 5 Fix fi < and denote by W^lO), W^{0) the stable and unstable manifolds 
of the fixed point U = of F. Consider the nonlinear map G = Fng o Fn^-i o ■ ■ ■ o 
F_nQ o ^-2no-i^ Equation ^9^) possesses an homoclinic orbit to if and only if 
W'{0) andGiWiO)) intersect. 

Lemma O is hard to use for analyzing homoclinic solutions since it involves a 
nonlinear transformation G instead of a linear one as in lemma |H However one can 
recover the linear case when replacing F„ by a suitable approximation F„, equal to 
Fn up to higher order terms. This is possible thanks to property ( I98l) of lemma [6] 
below. In the sequel we note 



DF(0) 



1 
-1 2-/i 
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Lemma 6 Consider the collection of maps (—tiq < n < uq) defined by 

F„, = A^FoA-\, (96) 
where A^no-i = I (^i^d for n > —no 

An = Ln Ln-l ■ ■ ■ L_no " L„ = (/ + Tn)L. (97) 

The map Fn is a leading order approximation of Fn, i.e. Fn = Fn + 0{6 \\{a,(3)\\^). 
Moreover one has the property 

Fn, O Fn,-1 O ... O F.no= A F^^°+\ (98) 

where A = An^ reads 

A = Ln, Ln,-i ■ ■ ■ L.no = / + 0{e). (99) 

Proof. 

First we note that the sequence An satisfies v4_„q = / + T_„q and 

An+l = iI + Tn+l)LAnL-' (100) 

for all n > —uq — 1. It follows for —uq < n < uq 

Fn = {I + Tn) L An-l L'' F O A'^. (101) 

Now let us note that An = I + 0{e). Moreover, the following identity holds true for 
any parameter-dependent matrix M G M2(M) with ||M|| = 0{e) 

F o{I + M) = L{I + M)L-^F + 0{e\\{a,p)f). (102) 

Consequently one has also 

F = L{I + M)L-^F o{I + M)-i + Oie \\ (a, f3) f ). 

Using this property in equation f llOip leads to 

Fn = {I + Tn)F + 0{e\\{a,P)f)- 

Using (lU^l) this yields F„ = F„ + 0(£: ||(a,/5)||^), therefore F„ is a leading order 
approximation of F„. Property (1981) follows directly from the definition of F„. 

□ 

It is worthwhile stressing that A = DG{0), where G is the nonlinear transfor- 
mation introduced in lemma [5l 

Now we fix in addition Fn = F = Fn for \n\ > + 1. According to lemma [6] 
we have also F„ = F„ + 0{e \\{a,P)\\^) for \n\ < uq. In the sequel we approximate 
system f l^ by the new one 

Un+l = FniUn). (103) 

Property (1981) implies the following result, since iy"(0) is invariant under i^2no+i_ 
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Lemma 7 Fix p < and denote by W^^O), ly" (0) the stable and unstable manifolds 
of the fixed point U = of F. Equation UU!^) possesses a solution Un homoclinic to 
if and only if W^*(0) and A{W^{Q)) intersect, where the matrix A = I + 0{e) is 
defined in lemma\^ The intersection point corresponds to f/no+i- 

Consequently, as in section 14.11 one recovers the problem of finding the intersec- 
tion of 14^'*(0) with the image of iy"(0) under the (near-identity) linear transforma- 
tion A. Note that A = J + Tq in the single defect case no=0. 

Here we shall not attempt to relate the bifurcations of breather solutions of ([6]) 
with the properties of the inhomogeneities, via an analysis of homoclinic solutions of 
( 11031) . This question will be considered in future works using the simplification pro- 
vided by lemma [71 As for a single defect, for i? < one can expect multiple tangent 
bifurcations between (deformations of) site-centered and bond-centered breathers as 
inhomogeneities are varied, due to the winding structure of W^"(0) and iy*(0). 

It is now interesting to compute the leading order contribution of the sequence 
of inhomogeneities to the matrix A. This is the object of the following lemma. 

Lemma 8 The matrix A of lemma takes the form A = I + M + 0{e'^ + e\fi\) , 
where 





n=0 



Ml 



12 




n=0 



21 




n=0 




n=0 



Pn 



fi2 



■n 



Proof. 



Since T„ = 0(e) it follows from definition (p9|) 




(104) 



n=0 
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Now we use the expansions 

T„ = M„ + 0(£2 + £|/i|), Mn ' ^ ° 

L = L, + 0{\^l\), L,= 



to obtain 



where 



1 
-1 2 

2no 



A = I + M + 0{e^ + e\fi\), M = ^ M„„_„ 4"", 



n=0 



in ^ f -n + l n 
^ \^ —n n + 1 

Then simple computations lead to the coefficients of M provided above. 

□ 

Interestingly, lemma [H] shows that the influence of the inhomogeneities on the 
set of homoclinic solutions depends at leading order (via the matrix / + M) on 
algebraically-weighted averages of {k} and {p}. 

Now let us return to the original parameters m„, dn, an, kn describing the lattice 
inhomogeneities (see equations (I5l) and (Ej)), with m„ = (i„ = a„ = for \n\ > uq + I, 
fc„ = for n < -no - 1 and n > no. Let us note e = ||{"^n/''^}||£oo + ll{'^n/'^}lkoo + 
ll{«n/o}||^oo + \\{kn/k}\\i^. One obtains 

h . — k 

k 

where 



dam 

is a linear combination of the on-site potential and mass defect impurities. Some 
coefficients of M can be simplified since 

2no 2no ^ rtQ-l 

^ ^ ^no-n = 0(5 ), ^ ^ Ti KnQ—n = ^ ^ kn -\- 0{£ ). 

n=0 n=0 n=— no 

Noting 

no "0-1 
n=— no n=— no 

one finally obtains A = I + M + O^e"^ + £\p\) with 
M = 



Mu -Mu + nolo - h 

Mil + {no + l)/o - h -Mil 
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and Mil = + l)/o - (St^o + 1)A + ^2 + -^o- 

Consequently, the matrix A depends (at leading order in e and fi) on the average 
values Jo, Jo of r„, kn/k, and on the weighted averages Ji, I2 of r„ (with linear and 
quadratic weights respectively). Since Tr(M) = 0, it follows Tt{A) = 2 + 0(e^+e|/i|) 
and Det(A) = l + 0(£^ + e|/i|). As a consequence, in order to study the spectrum of 
A for £, /X (and determine to which type of linear transformation it corresponds) 
it would be necessary to compute the quadratic terms in {e, fj,) in its expansion. 



5 Numerical results 

We have performed numerical computations in order to check the range of validity 
of the analysis of section 14.11 and in particular if discrepancies appear for large 
amplitude solutions or if parameters {mQ,uj) are moved away from (0,f2). More 
precisely, we have computed breather solutions of the Klein- Gordon lattice 

CJ^(1 + T^n)-^ + ^^V\yn) = Vn+l - 2?/„ + Un+l (105) 

with a single mass defect m„ = mo(5ri,o and periodic boundary conditions y-N(t) = 
yN(t)- In general we have used a lattice with 101 particles, except for the compu- 
tations of breathers with algebraic decay (case u = Q) where 401 particles have 
been considered. The computations have been compared with homoclinic orbits to 
of the two-dimensional map (17T1) . For the numerical computations we have always 
fixed f2 = 10 (recall Q is the lower phonon band edge for the infinite system). This 
can be done taking, for instance, k = 0.01, and d = 1 in the original problem 
For the potential V we have chosen a polynomial of degree 4 with V (0) = 1. 



5.1 Hard potentials 

To start we have considered the simplest case of a hard potential, i.e. a potential 
with a strictly positive hardening coefficient B (see definition ( 122|) ). We have chosen 

2 4 

V{x) = y + ^ , (106) 

for which B = 75. 

In this case, the reduced map ( 1711) possesses a unique orbit homoclinic to 
in the sector a > 0, (3 > 0, for ttiq > and uji < u < Q. An example of this 
homoclinic orbit is shown in figure for a frequency u = 9.9 (/i = —1.99) and a 
mass defect ttiq = 0.05. In figure [H] (left panel) we compare the approximate solution 
Vn = Pn COS t obtained with this homoclinic orbit (circles) with the exact breather 
profile computed with the standard numerical method based on the anti-continuous 
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limit [MA96j (continuous line). The agreement is excellent even if the solution profile 
is very localized. Indeed, as one computes the eigenvalues a, cr~^ (equation flTUl) ) of 
the linearized map ( I75l) with mo = 0, one obtains a ~ 0.27, which implies a strong 
spatial localization visible in figure [HI The accuracy of the centre manifold reduction 
(a priori expected for a ~ 1) is surprisingly good in this parameter regime. 




-10 -5 5 10 -10 -5 5 10 

n n 



Figure 8: Comparison between the profile of a breather solution (continuous line) 
of the Klein-Gordon system fllOSp with hard potential (11061) and the approximate 
solution Hn = PnCost (circles) constructed with the homoclinic orbit of ( ITTi) . We 
have considered a mass defect mo = 0.05. In the left panel we have chosen a 
frequency u = 9.9 (/x = —1.99). In the right panel we have fixed u = 9.837 
{n = —3.23) very close to ui (note the change of scale for the vertical axis). 

The breather solution can be continued for decreasing frequencies up to ui ~ 
9.8369, which is the frequency of the linear defect mode at which the breather 
solution bifurcates. Figure [H] (right panel) compares again the numerically computed 
breather profile and the approximate solution obtained with the homoclinic orbit, 
but now very close to this bifurcation point (at u = 9.837, i.e. fi = —3.23). We 
still observe an excellent agreement. Note that the oscillations amplitudes are very 
small, but the solution is still strongly localized. 

For increasing frequencies the continuation path ends up at the lower edge of 
the phonon band uj = Q (yU = 0). For this particular frequency value the breather 
solution (see continuous line in figure [HI left panel) presents an algebraic decay which 
is very well described by approximation ( 186|) . This approximation fails to describe 
the maximum amplitude of the oscillation /3o for these parameter values. This is not 
surprising since /?o is not small, and j3n varies rapidly near n = 0, hence mo should 
be further decreased to attain the domain of validity of the ansatz (1841) near the 
solution centre. 

However, the value of /3o obtained from the exact homoclinic orbit of ( |7T]) fits 
very well the maximum amplitude of the breather solution, as it is shown in figure [9l 
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right panel. Note that the agreement is very good even for very large amplitudes or 
very large mass defect i.e. far away from the parameter range in which the centre 
manifold reduction and normal form analysis are valid in principle. 




Figure 9: Left panel : breather solution at the lower edge of the phonon band on = 
{fi = 0) for a mass defect mo = 0.05 and the symmetric potential V{x) = x"^ /2+x^/4. 
The continuous line corresponds to the numerically computed breather solution. The 
circles represent approximation flSUj) of the homoclinic orbit that fits very well the 
algebraic decay of the breather tails. Right panel : the continuous lines represent 
the amplitude of the breather solution at n = (upper curve) and n = 1 (lower 
curve) versus mass defect. The circles correspond to the homoclinic solution of the 
nonlinear map ( 1711) (the upper plot represents (3o and the lower plot /?i). 



It is interesting to remark that the accuracy of this fit depends on the symmetry 
of the potential V{x) we have chosen. Figure fTOl shows what happens if we add to the 
polynomial potential f ll06p a cubic term x^/6 that breaks its symmetry. The range 
of validity of our leading order approximation reduces significantly. A similar result 
was obtained in reference |SRJCA04] for breather solutions in spatially homogeneous 
Fermi-Pasta-Ulam lattices. Obviously the agreement would be improved by taking 
into account the Taylor expansion of the reduction function (see theorem [3]) and 
computing the normal form at a higher order. 

Finally we have numerically checked that all breathers solutions in the gap ui < 
uj < Q are spectrally stable, at least for the value of the frequency parameter f2 = 10 
we have considered. This result is in agreement with the stability analysis performed 
in reference |CAR05j for Klein-Gordon chains with an impurity and at low coupling 
k. 
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Figure 10: Same computation as in figure [H right panel, but now for tlie asymmetric 
potential V{x) = + + x^/A. 



5.2 Soft potentials 

In the case of soft potentials (when the coefficient B defined by ( l22l) is strictly nega- 
tive), the situation is far more complex due to the much more intricate structure of 
the intersections between the stable and unstable manifolds. Therefore one expects 
a richer bifurcation scenario as parameters (breather frequency, mass defect) are 
varied. Our computations have been performed with the symmetric potential 



for which B = —75. 

Let us recall some basic features of the analysis performed in section 14.1.41 in 
order to compare the results with numerical computations. For the (truncated) 
reduced mapping (1711) with mo = 0, figure [7] shows some intersections of stable and 
unstable manifolds emanating from the saddle point at the origin, for a frequency 
value uj = 9.9 < Q. Iterating the map with an initial condition Ui at the homoclinic 
point with label 1, we obtain an homoclinic orbit which corresponds to a one-site 
breather centered at n = 0. With an initial condition f/i at the homoclinic point with 
label 2, the corresponding breather is a two-site breather with maximal amplitude 
at n = and n = 1. An initial condition Ui at the homoclinic point with label 3 
(symmetric of point 1 respect to the line a = (3) corresponds to a one-site breather 
centered at site n = 1. 

The dashed line of figure [7] depicts the image of the unstable manifold by the 
linear shear 74(a;,mo) for mo = 0.005. As mo increases A{u!,mo) W^{0) moves 
further down so that intersection points 2' and 3', corresponding to homoclinic orbits 
of the inhomogeneous problem, becomes closer and closer. So there exists a critical 



V{x) 



2 




(107) 
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value of mo for which these intersection points colhde and then disappear. In fact 
we have checked numerically that this tangent bifurcation occurs at a critical value 
mo G (0.00963,0.00964) for problem ( I7T1) . This critical value can be approximated 
using equations (l90l)-( l9Ti) . which yields uiq ^ 0.009632 in the present case. These 
results correspond very precisely to a breather bifurcation numerically observed in 
the Klein-Gordon chain (at a critical value tjiq G (0.00963,0.00964)) and depicted 
in figure [TTl 

The upper branch of figure [TTlfa) represents the energy of a breather solution 
corresponding to point 2'. For mo ~ 0, the breather has a maximal amplitude at 
sites 77, = 0, 1. A profile of this breather for mo = 0.0093, close to the bifurcation 
point, is shown in figure [TlT b). where the amplitude is now much larger at = 1. 
The lower branch of figure [TlTa) represents a one-site breather centered at n = 1 
and corresponds to point 3'. Its profile for mo = 0.0093 is shown in Figure fTlT c) . 
According to analytical results of reference |CAR05] . at low enough coupling k 
the solutions on the lower branch are spectrally stable whereas the solutions on 
the upper branch are unstable. We have checked this result numerically for the 
parameter values of figure [TT] by computing the Floquet spectra of the two families 
of breather solutions. 



As in section 15.11 we have also computed one-site breathers centered at the 
mass defect, corresponding to point 1' in figure [71 Again we have found an excellent 
agreement between the numerically computed breather profiles and the approximate 
solutions obtained using the map fl7T]) . As expected from the analysis of section H.1.4l 
these breathers survive up to mo = m;(u;), i.e. up to a much higher value of mo 
than the families 2', 3' described above. 

A part of the intersecting stable and unstable manifolds is shown in the left 
panel of figure [121 Due to their complicated windings, new intersections points 
appear between A^UjUIq) Vr"(0) and ly^(O) as mo is chosen in certain windows of 
the parameter space, giving rise to new homoclinic solutions of (1701) . 

An example is shown in the region marked with a rectangle (see the details 
in the right panel of figure [T2l) . For some value of mo G (0.01064,0.01065), a 
new intersection point between y4(ci;,mo) W^"(0) and iy**(0) appears. As mo is 
further increased, this inverse tangent bifurcation gives rise to two new homoclinic 
points 5' and 6'. Correspondingly, we have numerically checked that an inverse 
tangent bifurcation occurs in the Klein-Gordon chain at a critical value of mo G 
(0.01064,0.01065), giving rise to new breather solutions which do not exist in the 
homogeneous chain. 

The point 4' in figure [121 also exists for mo = 0. Returning to figure [71 it 
is obtained by applying the inverse map to the point with label 2. In the 
homogeneous limit mo = 0, this homoclinic point corresponds consequently to a 
two-site breather centered between n = 1 and n = 2. As figure [121 shows, an 
increase of the mass defect mo moves point 5' against point 4' until they collide 
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Figure 11: Tangent bifurcation between breather solutions numerically computed in 
a Klein-Gordon chain with a soft potential. The chain presents a mass defect mo at 
n = 0, and the bifurcation occurs as mo is increased. In the left panel, the breathers 
energies E = Zlnez ^^^(^/"(O)) + iVn+iiO) - y„(0))^/2 are depicted versus mo (the 
breathers are even in t with frequency uj = 9.9). For mo ~ 0, the upper branch 
represents a two-site breather centered between sites n = and n = 1. The lower 
branch represents a one-site breather centered at n = 1. The breathers profiles close 
to the bifurcation point are plotted in the right panels (the value of mo is marked 
with a dashed line in the left panel). 



and disappear through a new tangent bifurcation. This tangent bifurcation is also 
numerically found in the Klein Gordon chain at critical value of the mass defect very 
close to the theoretical one (in both cases one obtains mo ~ 0.01268). 

Figure [T3] shows the bifurcation diagram of the numerically computed breathers 
corresponding to homoclinic points 4', 5', 6' (left panel), and gives their profiles for 
a given value of mo in the right panel. A numerical computation of Floquet spectra 
shows that all these breathers are unstable. 

As a conclusion, we have seen that the truncated normal form (!70l) allows one 
to predict with a high precision certain breather bifurcations in the Klein-Gordon 
chain, which occur as the mass defect mo is varied. These bifurcations depend on 
the fine structure of the windings of the stable and unstable manifolds of the origin, 
computed on the truncated normal form without defect. 
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Figure 12: Emergence of new intersection points between A{uj, mo)W^ (dashed 
curve) and (drawn with a full line) as the mass defect is increased. The fig- 
ure corresponds to mo = 0.012 and u = 9.9. The right panel shows a zoom of the 
left panel over the region marked with a rectangle. The new homoclinic points 5' 
and 6' correspond to new breather solutions of the Klein-Gordon lattice. The point 
with label 4' corresponds to a two-site breather, which exists in the homogeneous 
lattice and persists for mo < 0.012. 



Upper branch 






A 












J. K. 










25 -20 -15 -10 


-5 5 

Central branch 


10 


15 


20 


25 


y , V 


25 -20 -15 -10 


-5 5 

Lower branch 


10 


15 


20 


25 
















"v. 











-25 -20 -15 -10 -5 5 10 15 20 25 



Figure 13: Bifurcation diagram of breather solutions numerically computed in the 
Klein-Gordon chain, with a soft potential and a mass defect mo at n = 0. In the left 
panel the breathers energies E are depicted versus mo (see the definition of E in the 
caption of figure fTTl) . The breathers frequency is a; = 9.9. The lower branch at the 
left of the vertical line corresponds to a two-site breather centered between n = 1 
and n = 2. The right panel shows the profiles of the three breathers for mo = 0.012, 
when all of them coexist (the value of mo is marked with a vertical line in the left 
panel) . 
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(UMR 5219) where a part of this work has been carried out during a visit in Sept.- 
Oct. 2006. 
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